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ABSTRACT 


A  high-resolution,  multi-level,  primitive  equation  ocean  model  is  used  to 
examine  the  response  of  an  idealized,  flat-bottom,  eastern  boundary  oceanic 
regime  on  a  beta-plane  to  climatological  average  (1980-1989),  individual  year, 
and  multiple  year  wind  forcing.  The  focus  of  this  study  is  the  California  Current 
System  along  the  coastal  region,  from  35°  N  to  47.5°  N,  off  the  West  Coast  of 
North  America. 

Two  types  of  experiments  are  conducted.  The  first  type  forces  the  model 
from  rest  with  climatological,  1981,  and  1983  monthly  winds  to  examine  the 
generation  phase  of  features  such  as  currents,  upwelling,  meanders,  eddies,  and 
filaments.  The  second  type  continues  the  forcing  from  the  previous  years  to 
examine  the  maintenance  of  these  features. 


In  the  first  type  of  experiments,  the  following  features  are  observed:  a 
poleward  coastal  surface  current  near  the  start  and  end  of  each  year,  an 
equatorward  surface  current,  a  poleward  undercurrent,  upwelling,  meanders,  and 
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I .  INTRODUCTION 


The  California  Current  System  (CCS)  is  an  eastern 
boundary  current  (BBC)  off  the  west  coast  of  the  United 
States  of  America.  A  primitive  equation  numerical  model 
developed  by  Batteen  gjfc.  al.  (1989)  is  used  to  simulate  the 
currents,  eddies  and  filaments  in  the  CCS.  In  this  study 
monthly  winds  from  the  Buropean  Centre  for  Medium- Range 
Weather  Forecasts  (ECMWF)  global  data  set  from  1980  to  1989 
(which  include  El  Nifio  events)  are  used  to  force  the  model. 

The  thesis  is  organized  as  follows:  Chapter  II 
describes  the  region  modeled,  the  El  Nifio -Southern 
Oscillation  (ENSO) ,  and  some  previous  modeling  work  in  this 
region.  The  numerical  model  which  was  used  here  along  with 
the  experimental  setup  are  presented  in  Chapter  III.  The 
results  of  the  various  model  simulations  are  presented  in 
Chapter  IV.  Chapter  V  summarizes  the  work,  discusses  naval 
relevance,  and  presents  recommendations  for  future  research. 
Comparisons  of  model  simulations  with  available  observations 
are  shown  in  both  Chapters  IV  and  V.  Six  appendices  follow 
which  contain  the  detailed  results  of  the  barotropic  and 
baroclinic  instability  analyses  for  the  main  experiments . 
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II .  BACKGROUND 


A.  BRIEF  DESCRIPTION  OF  REGION 

The  west  coast  of  North  America  serves  as  an  eastern 
boundary  for  the  circulation  of  water  in  the  North  Pacific 
Ocean.  The  entire  flow  along  the  coast  is  referred  to  as 
the  CCS,  which  is  a  classical  subtropical  EBC  system. 

The  CCS  consists  of  three  separate  component  flows,  each 
of  which  varies  both  spatially  and  seasonally  (see  Figure 
2.1):  an  equatorward  surface  current,  a  poleward 
undercurrent  and  a  poleward  surface  current  adjacent  to  the 
coast  (Tisch  fit  al. .  1992).  This  latter  current  can  be  a 
surface  manifestation  of  the  undercurrent  or  a  separate 
current . 

The  predominant  flow  is  am  equatorward  surface  flow 
called  the  California  Current  (CC) .  The  second  basic  flow 
is  the  poleward  undercurrent  which  has  been  called  both  the 
California  Undercurrent  (CUC)  (e.g.,  by  Batteen  fit  al. . 

1989)  amd  the  California  Counter  Current  (CCC)  (e.g.,  by 
Hickey,  1979) .  This  is  a  relatively  weak  subsurface  flow 
(around  5  cm  s'1) . 

The  third  basic  flow  is  the  Inshore  Current  (IC) ,  which 
is  known  as  the  Davidson  Current  (DC)  north  of  Pt. 
Conception.  This  is  also  a  relatively  weaik  poleward  flow 
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found  at  the  surface  and  near  the  coast.  There  are  other 
ICs  known  as  the  Southern  California  Countercurrent  (SCC)  to 
the  south,  and  the  Southern  California  Eddy  (SCE)  inshore  of 
the  Channel  Islands  within  the  California  Bight  (Hickey, 
1979),  which  are  both  equatorward  of  the  model  region  (see 
below)  to  be  studied. 

The  CC  is  a  broad,  slow  (~10  cm  s'1)  flow  which  can 
extend  offshore  around  1000  km  (Hickey,  1979) .  Superimposed 
on  this  mean  flow  are  highly  energetic  mesoscale  eddies  and 
meandering  jets  which  can  have  peak  velocities  near  80  cm 
s"1  (Mooers  and  Robinson,  1984).  Figure  2.2  depicts  these 
features . 

A  section  of  the  CCS  from  35°  to  47.5°  N  and  extending 
12.5  degrees  offshore  from  the  coast  will  be  the  model  study 
area  and  is  depicted  in  Figure  2.3.  Major  observational 
programs  which  have  sampled  this  region  are  quite  numerous 
and  include:  California  Cooperative  Oceanic  Fisheries 
Investigation  (CalCOFI) ,  Coastal  Ocean  Dynamics  Experiment 
(CODE)  I  and  II,  Large-Scale  West  Coast  Shelf  Experiment 
(nicknamed  SuperCODE) ,  Coastal  Upwelling  Experiment  (CUE)  I 
and  II,  Ocean  Prediction  Through  Observation,  Modeling  and 
Analysis  program  (OPTOMA) ,  and  Central  California  Coastal 
Circulation  Study  (CCCS) .  Figure  2.4  depicts  some  of  the 
areas  of  observations. 
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B.  CLIMATOLOGICAL  WINDS 


The  North  Pacific  Subtropical  High  dominates  the 
atmospheric  circulation  off  the  west  coast  of  the  United 
States  of  America,  and  varies  in  position  and  strength  by 
season  (Huyer,  1983) .  In  the  winter  (January)  the  high  is 
about  1020  mb  and  centered  near  30°  N,  135°  W  (Figure  2.5a). 
In  the  summer  (July)  the  pressure  increases  to  1025  mb  and 
the  center  moves  northward  to  38°  N,  150°  W  (Figure  2.5b). 
Over  the  continental  west  coast  the  typical  pressure  is  1020 
mb  in  the  winter,  and  1015  mb  during  the  summer,  the  latter 
in  response  to  a  thermal  low  (of  1005  mb)  in  the  central 
valley  of  California  (Huyer,  1983) .  This  results  in  a 
tighter  pressure  gradient  against  the  west  coast  during  the 
summer.  As  this  gradient  tightens,  stronger  winds  can 
result. 

Propagating  atmospheric  disturbances  (Halliwell  and 
Allen,  1987) ,  along  with  other  atmospheric  mesoscale 
phenomena  (Huyer,  1983) ,  further  complicate  the  wind  pattern 
on  a  daily,  as  well  as  on  a  seasonal  basis  along  the  coast. 
The  winds  are  also  affected  by  coastal  atmospheric  boundary 
layer  processes  within  100  to  200  km  of  the  coast  (Halliwell 
and  Allen,  1987) . 

The  oceanic  response  to  summer  and  winter  wind  patterns 
are  most  notable  near  the  coast.  The  equatorward  winds  of 
summer  can  produce  nearshore  upwelling,  while  the  divergent 
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wind  flow  of  winter  can  force  downwelling  in  sane  areas 
(Rienecker  and  Mooers,  1986) . 

The  shifting  of  the  Pacific  High,  along  with  the 
Aleutian  Low  moving  to  the  southeast,  can  cause  the  winds  to 
diverge  near  Cape  Mendocino  {located  at  about  40°  N)  during 
the  winter.  North  of  the  Cape,  northerly  winds  are 
encountered  as  the  wind  blows  around  the  eastern  side  of  the 
Aleutian  Low,  while  south  of  the  Cape,  southerly  winds 
result  as  wind  flows  around  the  North  Pacific  High. 

The  summer  wind  pattern  is  established  by  the 
interaction  between  the  North  Pacific  High  and  the  thermal 
low  in  the  southwest  United  States  (Nelson,  1977;  Halliwell 
and  Allen,  1987) .  Winds  are  predominately  equatorward  on 
both  sides  of  Cape  Mendocino  during  the  summer  (Nelson, 

1977) . 

Climatological  wind  stress  fields,  shown  in  Figure  2.6 
for  the  months  of  January,  April,  July,  and  October,  depict 
regions  of  maximum  wind  stress  in  the  CCS  as  shaded  areas. 
The  strongest  wind  stress  is  located  off  Point  Conception  in 
April  and  migrates  northward  until  it  is  at  a  maximum  in 
July  off  Cape  Mendocino.  The  average  summer  wind  stress 
close  to  the  shore  is  equatorward  and  is  thus  favorable  for 
upwelling  at  the  coast  (Halliwell  and  Allen,  1987;  Nelson, 
1977;  Wickham  ££  al. .  1987). 
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C.  EL  Nlffo- SOUTHERN  OSCILLATION 

One  large-scale  two  way  coupling  of  the  atmosphere  and 
ocean  has  been  studied  for  many  years  and  is  referred  to  as 
the  El  Nifio- Southern  Oscillation  (ENSO) .  El  Nifio  refers  to 
the  wanning  of  Eastern  Pacific  ocean  waters,  while  the 
Southern  Oscillation  is  an  atmospheric  pressure  fluctuation 
between  the  Indian  Ocean  and  the  Eastern  Tropical  Pacific 
Ocean  (Philander,  1990).  Figure  2.7  illustrates  key 
properties  of  the  teleconnection  including  upper  troposphere 
and  surface  pressure  anomalies.  Note  the  shift  in  the 
Eastern  Pacific  Subtropical  High  during  an  ENSO  event. 

Explanations  for  this  phenomena  involve  both  an  oceanic 
and  atmospheric  teleconnection.  Rienecker  and  Mooers  (1986) 
explain  that  the  oceanic  connection  occurs  when  eastward 
propagating  equatorial  Kelvin  waves  turn  poleward  upon 
impacting  the  west  coast  of  the  Americas.  In  contrast,  the 
atmospheric  teleconnection  links  the  tropical  sea  surface 
temperature  (SST)  anomalies  to  the  Aleutian  Low  by  momentum 
transfer  through  an  intensified  Hadley  Cell  circulation 
(Bjerknes,  1966;  Philander,  1990) .  A  stronger  poleward  wind 
stress  subsequently  develops  at  the  coast  near  36°  N  when 
the  Aleutian  Low  deepens  and  moves  to  the  southwest. 
Simultaneously  the  North  Pacific  High  weakens  and  moves 
offshore. 

The  El  Niflo  of  1982  -  1983  was  one  of  the  strongest  on 
record.  The  Aleutian  Low  was  much  deeper  than  normal  during 


7 


the  winter  and  wind  stress  anomalies  were  at  a  maximum  in 
the  region  between  35°  N  and  45°  N  (Trenberth  fit,  al . .  1990). 
Anomalous  conditions  were  observed  from  San  Diego  to 
Vancouver  Island  (Rienecker  and  Mooers,  1986) .  There  were 
two  anomalous  peaks  in  the  SST,  the  first  occurring  in 
February  to  March  of  1983  off  northern  California,  while  the 
second  one  was  evident  in  September  to  October  of  1983  near 
Point  Conception  (Rienecker  and  Mooers,  1986) ,  the  latter 
equatorward  of  our  model  domain. 

D.  NUMERICAL  MODEL  STUDIES 

There  have  been  numerous  models  used  over  the  past  few 
decades  to  model  EBCs.  particularly  the  CCS.  Extensive 
reviews  can  be  found  in  Allen  (1980) ,  Chelton  (1984) ,  and 
O'Brien  fit  al*  (1977) .  No  attempt  is  made  to  provide  a 
comprehensive  review  of  the  models  here;  instead,  a  summary 
of  modeling  work  that  has  led  directly  to  this  study  is 
presented . 

Early  work  included  that  of  Pedlosky  (1974)  and 
Philander  and  Yoon  (1982)  who  used  steady  wind  stress  and 
transient  wind  forcing,  respectively.  Carton  (1984)  and 
Carton  and  Philander  (1984)  investigated  the  response  of 
reduced  gravity  models  to  realistic  coastal  winds.  At 
periods  exceeding  100  days  the  alongshore  currents  began  to 
weaken  and  disperse  as  a  series  of  alternating  jets. 
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Another  series  of  experiments  were  done  by  McCreary  et 
al .  (1987)  utilizing  a  linear  model  with  both  transient  and 
steady  wind  forcing  with  Laplacian  diffusion  in  which  no 
eddies  or  filaments  developed.  This  work  was  extended  using 
a  full  primitive  equation  model  by  Tielking  (1988)  who 
obtained  similar  results. 

Batteen  fit  al.  (1989)  used  a  primitive  equation,  multi¬ 
level  model  with  biharmonic  heat  and  momentum  diffusion 
which  had  both  steady  forcing  constant  throughout  the  domain 
and  meridionally  varying  forcing.  In  both  cases,  when  the 
baroclinic  shear  between  the  surface  current  and 
undercurrent  became  strong  enough,  eddies  and  filaments 
developed.  This  was  the  first  BBC  model  to  simulate 
currents,  eddies  and  jets  in  the  CCS. 

The  work  of  Batteen  fit  al.  (1989)  was  extended  by 
Mitchell  (1993)  and  Haines  (1994) .  Mitchell  (1993)  forced 
the  model  with  different  sets  of  averaged  winds  and  obtained 
reasonable  results  utilizing  a  set  of  seasonally  varying 
full  climatological  wind  fields.  Haines  (1994)  enlarged  the 
model  domain,  improved  resolution  to  8  x  13  km  and  used 
interannual  vice  steady  or  seasonal  winds. 

This  effort  builds  on  the  model  used  by  Batteen  fit  al. 
(1989)  and  the  work  of  Mitchell  (1993)  and  Haines  (1994) . 

The  computer  code  for  the  model  has  been  modified  to  include 
a  barotropic  component,  while  the  model  resolution  has  been 
refined  to  8  x  ll  km  grid  spacing.  Seasonal  and  interannual 
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wind  forcing  was  used  to  force  the  model.  Model  simulation 
times  have  been  expanded  from  one  year,  as  in  Mitchell 
(1993)  and  Haines  (1994) ,  to  two  to  four  years  for  each 
model  run.  The  enhanced  simulation  times  should  allow 
analysis  of  the  maintenance  of  currents  and  eddies  as  well 
as  the  generation  of  these  features. 
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Longitude  (degrees  west) 


Figure  2.1  Generalized  circulation  schematic  of  the 
California  Current  System  (CCS) .  The  broad,  slow  surface 
equatorward  California  Current  (CC)  overlies  the  poleward 
California  Undercurrent  (CUC)  along  with  the  Inshore  Current 
(IC) ,  also  known  as  the  Davidson  Current  (DC) ,  Southern 
California  Countercurrent  (SCC) ,  and  Southern  California 
Eddy  (SCE) . 
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Figure  2.2  Section  of  California  Current  System  with 
schematized  positions  of  jets  and  synoptic- scale  eddies;  C, 
cyclonic  eddy;  A,  anticyclonic  eddy.  Arrows  indicate 
principal  surface  flows;  solid  arrows  are  definite  and  open 
arrows  are  inferential.  Also  shown  is  a  qualitative 
analysis  of  satellite  infrared  image  gray  shades  for  sea- 
surface  temperature  (SST)  predominant  patterns  on  l  August 
1982:  coolest  SST,  crosses;  warmest  SST,  dots  (from  Mooers 
and  Robinson,  1984) . 
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Figure  2.3  Domain  of  model  off  the  western  coast  of  the 
United  States  of  America.  Domain  bounded  by  35®  to  47.5°  N 
and  from  the  coast  to  12.5  degrees  offshore.  Model  covers 
an  area  1024  km  (1408  km)  in  the  cross -shore  (alongshore) 
direction  and  the  cross -shore  (alongshore)  horizontal 
resolution  is  8  km  (11  km) . 
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Figure  2.4  Map  of  a  portion  of  the  west  coast  of  North 
America,  showing  the  location  of  some  of  the  past  studies  in 
the  CCS  (from  Buyer  fit  al*  #  1989). 
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Figure  2.5  Long-term  mean  atmospheric  pressure  at  sea  level 
for  (a)  January  and  (b)  July  (from  Huyer,  1983) . 
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2.6  Long-term  (1854-1972)  mean  wind  stress  (dynes 
cm’  )  by  one -degree  squares,  from  ship  reports  (Nelson, 
1977).  Contours  correspond  to  constant  magnitudes  of  0.5, 
1.0,  and  1.5  dynes  cm'2;  areas  where  the  wind  stress  exceeds 
1.0  dynes  cm'2  are  shaded  (from  Huyer,  1983). 
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Figure  2.7  An  illustration  of  the  distributions  of  key 
properties  involved  in  the  teleconnection  process  during  a 
typical  El  Nilio- Southern  Oscillation  event.  (Left)  Upper 
troposphere  anomalies  of  geopotential  height  for  the 
Northern  Hemisphere  winter.  (Right)  Contours  of  surface 
barometric  pressure,  with  surface  wind  anomalies  (arrows)  in 
the  tropics  (from  Enfield,  1989) . 
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III.  MODEL  DESCRIPTION  AND  EXPERIMENTAL  DESIGN 


A.  MODEL  DESCRIPTION 

1.  Basic  Model 

The  model  used  in  this  research  is  a  slightly 
modified  version  of  that  used  by  Batteen  fit  &1-  (1989,  1992) 
which  was  adapted  from  Haney's  (1974,  1985)  closed  basin 
model.  This  limited  area  BBC  model  is  multi-level,  uses 
non-adiabatic  primitive  equations  on  a  0  plane,  and  has  both 
barotropic  and  baroclinic  components .  The  following 
approximations  are  made:  hydrostatic,  Boussinesq,  and  rigid 
lid.  (For  a  more  detailed  summary  of  the  model  and  the 
governing  equations,  refer  to  Batteen  fit  al.  (1989)). 

2 .  Friction 

The  type  of  diffusion  used  can  be  important  for  the 
generation  of  mesoscale  eddies.  Biharmonic  rather  them 
Laplacian  lateral  heat  and  momentum  diffusion  is  used  in 
this  study.  The  use  of  biharmonic  lateral  diffusion  allows 
mesoscale  eddy  generation  via  barotropic  (horizontal -shear) 
or  baroclinic  (vertical -shear)  instability  mechanisms,  or 
both,  as  Holland  (1978)  and  Holland  and  Batteen  (1986)  have 
shown.  The  biharmonic  coefficients  are  the  same  (i.e.,  2.0 
X  1017  ernes'1)  as  used  in  Batteen  fit  fil.  (1989).  Weak  (0.5 
cm2s'1)  vertical  eddy  viscosities  and  conductivity  are  used. 
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A  simplified  quadratic  drag  law  is  used  to  parameterize 
bottom  stress  (Weatherly,  1972) . 

3 .  Boundary  Conditions 

The  domain  of  the  model  on  the  eastern  boundary  is 
closed  and  is  modeled  as  a  straight,  vertical  wall.  A  no¬ 
slip  condition  on  the  alongshore  velocity  is  used.  The 
kinematic  boundary  condition  of  no  flow  through  the  eastern 
boundary  is  imposed  on  the  cross -shore  velocity  component. 

A  modified  version  of  the  radiation  boundary  conditions  of 
Camerlengo  and  O'Brien  (1980)  is  used  for  the  open 
boundaries  to  the  south,  west,  and  north. 

At  the  surface,  the  rigid  lid  approximation  filters 
out  external  gravity  waves  and  allows  a  longer  time  step. 

On  the  ocean  bottom,  no  vertical  heat  flux  or  vertical 
velocity  is  allowed. 

4.  Domain  and  Resolution 

The  domain  of  the  model  is  an  elbow  shaped  region 
off  California's  coast  from  35°  N  to  42.5°  N  as  shown  in 
Figure  2.3.  This  covers  an  area  1408  km  alongshore  and  1024 
km  in  the  cross -shore  direction.  The  horizontal  resolution 
is  8  km  in  the  cross -shore  direction  and  11  km  in  the 
alongshore  direction.  Mesoscale  features  within  the  CCS 
should  be  resolvable  with  this  resolution  because  typical 
eddy  wavelengths  are  around  100  km  (Breaker  and  Mooers, 

1986) .  Bottom  topography  and  variations  in  the  coastline 
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have  been  omitted  to  focus  on  the  role  played  by  wind 
forcing  in  the  CCS. 

5.  Finite  Difference  Scheme 

A  staggered  grid  model,  using  the  Arakawa  and  Lamb 
(1977)  B-scheme  is  used  in  the  horizontal.  Variables  are 
defined  as  follows:  u,  v,  and  w  are  velocity  components 
eastward,  northward  and  vertical,  respectively,  t  denotes 
time,  T  is  temperature,  p  is  density,  p  is  pressure,  and  ^ 
is  the  streamf unction.  The  center  of  each  grid  point 
defines  u  and  v  while  T,  p,  w,  and  p  are  defined  at  the 
corners.  The  B-scheme  is  used  because  it  gives  a  superior 
representation  of  the  geostrophic  adjustment  process 
(Batteen  and  Han,  1981) . 

The  model  uses  ten  layers  in  the  vertical,  with 
constant  z-levels  at  depths  of:  13,  46,  98,  182,  316,  529, 
870,  1416,  2283,  and  3656  m.  This  spacing  scheme 
concentrates  more  layers  above  the  thermocline  in  the 
dynamically  active  portion  of  the  ocean,  consistent  with 
Haney  (1974) .  Variables  are  staggered  in  the  vertical  with 
u,  v,  p,  and  T  defined  at  each  level,  while  w  and  p  are 
defined  between  the  levels  which  lie  between  the  ocean 
surface  (w*0)  and  at  the  bottom  (z--H,  where  H  is  the 
depth) .  A  Matsuno  time  step  followed  by  ten  leapfrog  time 
steps  is  used  and  continuously  repeated  during  the  model 
run. 


21 


1 


6.  Method  of  Solution 

The  method  of  solution  is  straightforward  with  the 
rigid  lid  and  flat  bottom  assumption  because  the  vertically 
integrated  horizontal  velocity  is  subsequently  non- 
divergent.  The  vertical  mean  flow  can  be  described  by  a 
streamfunction  ^  which  can  be  predicted  from  the  vorticity 
equation,  while  the  vertical  shear  currents  can  be  predicted 
after  the  vertical  mean  flow  is  subtracted  from  the  original 
equations.  The  other  variables  T,  p,  w,  and  p  can  be 
explicitly  obtained  from  the  thermodynamic  energy  equation, 
equation  of  state,  continuity  equation,  and  hydrostatic 
relation,  respectively. 

7 .  Initial  Conditions 

An  exponential  temperature  profile  with  a  vertical 
length  scale  of  h-450  m  is  used  as  the  initial  mean 
stratification.  The  bottom  temperature  is  assumed  to  be  2° 

C  and  the  change  in  temperature  between  the  bottom  and 
surface  is  13°  C.  This  profile  was  derived  by  Blumberg  and 
Mellor  (1987)  and  has  been  used  by  Batteen  (1989)  and 
Batteen  ££  al.  (1989) . 

B.  WIND  AND  ATMOSPHERIC  FORCING 

The  winds  used  to  force  the  model  are  from  the  European 
Centre  for  Medium- Range  Weather  Forecasts  (ECMWF)  (Trenberth 
et  al . ,  1990)  and  have  been  provided  by  the  National  Center 
for  Atmospheric  Research  (NCAR) .  Winds  were  extracted  from 
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world  wide  mean  monthly  wind  velocities  that  were  provided 
at  2.5  degree  spacing  for  the  years  1980  to  1989  along  with 
120  month  climatology  for  the  area  of  study.  This  data  was 
linearly  interpolated  in  space  and  then  in  time  to  provide 
daily  forcing  at  the  desired  grid  points. 

The  total  heat  flux  across  the  surface  is  initially  set 
to  zero  by  choosing  an  initial  air  temperature  that  forces 
the  net  flux  of  longwave  radiation,  sensible  heat,  and 
latent  heat  to  zero,  as  in  Batteen  fit  &1-  (1989) .  The  air 
temperature  is  then  used  in  the  model.  Any  subsequent 
surface  heat  flux  forcing  is  thus  a  secondary  effect  of  the 
changes  to  SST  due  to  the  wind  and/ or  thermal  structure  of 
the  ocean. 

C.  SPECIFIC  EXPERIMENTAL  DESIGN 

Climatological  winds  were  used  for  365  days  (Experiment 
1)  and  then  used  to  restart  the  model  to  run  climatological 
winds  for  another  365  days  (Experiment  4) .  Next,  1981  winds 
were  run  for  365  days  (Experiment  2)  and  used  to  initialize 
the  run  of  365  days  with  1982  winds  (Experiment  5) .  Then 
1983  winds  from  a  cold  restart  were  run  for  365  days 
(Experiment  3)  and  then  used  to  initialize  the  run  of  366 
days  (leap  year)  with  1984  winds  (Experiment  6) .  Two 
additional  extended  runs  were  conducted  in  which  the  model 
ran  for  four  years.  Experiment  7  used  the  climatological 
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winds  four  times,  while  Experiment  6  used  the  monthly  winds 
of  1981  through  1984. 

D.  ENERGY  ANALYSIS  TECHNIQUES 

The  energy  technique  used  is  the  same  as  that  used  and 
described  in  Batteen  et  al-  (1992),  and  is  based  on  that  of 
Han  (1975)  and  Semtner  and  Mintz  (1977) .  This  is  done  to 
gain  a  better  understanding  of  the  flow  within  the  CCS  and 
of  the  types  of  energy  transfer  during  unstable  flow.  A 
brief  summary  of  the  method  follows. 

Kinetic  energy  (K)  is  calculated  for  the  horizontal 
components .  After  quasi -steady  state  is  reached  where  the 
total  kinetic  energy  is  nearly  constant,  mean  and  eddy 
kinetic  energies  are  calculated  using  the  averaged  stun  of 
squared  mean  and  eddy  horizontal  fields,  respectively. 

Next,  the  available  potential  energy  is  calculated  and  used 
to  determine  when  a  quasi -steady  state  is  reached  and  when 
statistics  should  be  collected.  Then,  both  mean  and  eddy 
available  potential  energy  are  computed. 

The  energy  transfers  follow  that  of  Semtner  and  Mintz 
(1977) ,  and  are  used  to  argue  for  the  type  of  instability 
mechanism  (e.g.,  barotropic,  baroclinic,  or  mixed)  which 
leads  to  the  initial  eddy  generation  in  each  experiment. 
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IV.  RESULTS  OF  EXPERIMENTS 


Experiments  1  through  3  study  the  effects  of  forcing  a 
model  ocean  domain  that  is  initially  at  rest.  Experiment  1 
uses  monthly  climatological  average  (1980  -  1989)  wind 
fields,  while  Experiments  2  and  3  use  monthly  average  wind 
fields  for  1981  and  1983,  respectively.  Experiments  4 
through  6  study  the  effects  of  forcing  a  model  domain  that 
has  been  spun  up  by  the  wind  field  from  the  previous  year. 
Experiment  4  uses  the  results  of  Experiment  1  as  the  start 
for  a  second  year  of  climatological  average  wind  fields, 
while  Experiment  5  and  6  use  the  results  of  Experiment  2  and 
3  as  the  start  for  1982  and  1984,  respectively.  Experiments 
7  and  8  use  the  results  of  longer  (4 -year)  time  runs  from 
climatological  and  1981  -  1984  wind  fields,  respectively,  to 
examine  the  maintenance  of  features  after  initial  spin  up. 

A.  EXPERIMENTS  STARTED  FROM  REST 

1.  Experiment  1  -  Climatological  Minds 

Experiment  l  was  forced  with  monthly  averaged 
climatological  (1980  -  1989)  winds.  Figure  4.1  depicts  the 
monthly  wind  patterns  over  the  model  domain  starting  with 
January  15  (day  15)  and  ending  on  December  15  (day  345) . 

The  atmospheric  pressure  pattern  for  January  has  a  low  to 
the  north  and  a  high  to  the  south  which  results  in  a  wind 
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divergence  near  40°  N  (Figure  4.1a).  This  pattern  of 
poleward  winds  north  of  40°  N  and  equatorward  winds  to  the 
south,  continues  through  February  and  March  (Figures  4.1b 
and  c) .  During  April  and  May  (Figures  4. Id  and  e) ,  the 
divergence  in  the  wind  flow  is  observed  to  migrate  poleward. 
By  June  (Figure  4. If),  an  equatorward  component  in  the  wind 
field  is  observed  along  the  entire  model  domain.  July, 
August,  and  September  (Figures  4.1g  through  i)  depict  the 
strongest  and  most  equatorward  winds.  By  October  (Figure 
4 . l j ) ,  the  winds  start  to  weaken  throughout  the  domain  and  a 
divergence  in  the  wind  flow  is  observed  in  the  north.  This 
divergent  wind  flow  pattern  continues  through  November 
(Figure  4.1k).  By  December  (Figure  4.11)  the  wind 
divergence  has  returned  to  around  40°  N. 

In  response  to  the  prevailing  poleward  winds  in  the 
north,  a  relatively  weak  poleward  coastal  surface  current  of 
~3  cm  s"1  was  observed  by  model  day  30  (Figure  4.2)  .  It 
retreated  further  northward  until  by  day  60  (not  shown)  it 
was  not  discernible.  An  equatorward  surface  current  was 
observed  along  most  of  the  coast  by  model  day  90  (Figure 
4.3)  due  to  the  prevailing  equatorward  winds  over  most  of 
the  model  domain.  A  poleward  undercurrent ,  caused  by  the 
meridional  variation  of  the  wind  field  (Batteen  et  al. . 
1989),  first  became  evident  at  model  day  60  (Figure  4.4), 
and  was  fully  established  by  model  day  90  (Figure  4.5). 
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Along  with  the  equatorward  surface  current, 
upwelling  also  occurred  and  was  discernible  by  day  120 
(Figure  4.6).  Meanders  in  the  equatorward  surface  current 
formed  by  model  day  210  and  are  evident  in  both  the  velocity 
and  temperature  contour  plots  (Figures  4.7  and  4.8, 
respectively) .  By  model  day  240,  meanders  began  to  form 
cold  core,  cyclonic  eddies  as  seen  in  the  temperature 
contour  and  velocity  plots  (Figures  4.9  and  4.10, 
respectively)  due  to  vertical  and/or  horizontal  shear 
instabilities  between  the  equatorward  jet  and  the  poleward 
undercurrent.  The  eddies  initially  appear  to  be  about  30  km 
(Rossby  radius  of  deformation)  in  diameter  and  extend  about 
100  km  off  the  coast,  and  appear  during  the  period  of 
maximum  wind  stress  (July  -  September)  . 

The  coldest  temperatures  due  to  upwelling  were 
reached  by  model  day  240.  After  model  day  270  no  upwelling 
was  discernible  (not  shown) .  At  the  end  of  the  year, 
meanders  and  eddies  in  a  weak  equatorward  surface  current 
were  still  present  (Figures  4.11  and  4.12).  A  very  weak 
poleward  coastal  surface  current  (not  shown) ,  and  a  poleward 
undercurrent  were  also  present  at  the  end  of  the  year 
(Figure  4.13) .  Note  that  no  filaments  developed  during  this 
first  climatological  year. 

Except  for  the  lack  of  development  of  filaments, 
these  model  results  are  consistent  with  available 
observations  cited  in  the  literature.  For  example,  the 
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eddies  have  rotational  velocities  of  around  50-80  cm  s*1 
which  is  in  agreement  with  observations  of  Mooers  and 
Robinson  (1984) .  The  poleward  undercurrent  is  relatively 
weak  (~5  cm  s*1  or  less)  and  is  in  agreement  with  available 
observations  by  Huyer  fit  &1-  (1989) . 

A  time  series  of  total  kinetic  energy  per  unit  mass 
over  the  entire  domain  is  shown  in  Figure  4.14.  In  this 
figure  there  are  around  three  periods  when  the  kinetic 
energy  becomes  quasi -steady.  These  time  frames  correspond 
to  model  days:  207  -  216,  228  -  231,  and  297  -  303.  Since 
meanders  and/or  eddies  appeared  during  these  time  periods, 
more  detailed  analysis  was  performed  on  these  intervals  to 
determine  the  type  of  instability  mechanism  that  could 
generate  these  features.  Barotropic  instability  can  result 
from  horizontal  shiar  while  baroclinic  instability  can 
result  from  vertical  snear  in  the  current.  As  Figure  4.12 
shows,  there  is  considerable  horizontal  and  vertical  shear 
between  the  equatorward  surface  current  and  the  poleward 
undercurrent.  As  a  result,  both  types  of  instabilities 
(mixed) ,  can  be  present  simultaneously.  Energy  transfer 
calculations  which  consist  of  barotropic  (mean  kinetic 
energy  to  eddy  kinetic  energy)  and  baroclinic  (eddy 
potential  energy  to  eddy  kinetic  energy)  components  were 
performed  on  the  three  time  frames  above. 

The  results  for  the  instability  analyses  for  these 
periods  are  summarized  in  Table  4.1  (see  Appendix  A  for  the 
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energy  transfer  plots  for  this  experiment) ,  and  show  the 
following:  For  model  days  207  to  216,  228  to  231,  and  297 
to  303  both  barotropic  and  baroclinic  (i.e.,  mixed) 
instabilities  were  present  in  the  coastal,  equatorward 
region  of  the  domain.  Additional  analysis  of  model  days  297 
to  303  shows  barotropic  instabilities  to  be  weak  (Figure 
A. 5)  and  baroclinic  instabilities  to  be  dominant  (Figure 
A. 6)  during  this  period. 

2.  Experiment  2  -  1981  Winds 

Experiment  2  was  forced  with  the  monthly  ECMWF  winds 
from  1981.  Figure  4.15  depicts  the  monthly  wind  patterns 
within  the  model  domain  starting  with  January  and  running 
through  December.  Unlike  the  climatological  winds,  January 
starts  off  with  poleward  flow  throughout  the  model  domain 
(Figure  4.15a).  February  (Figure  4.15b)  shows  the  typical 
wind  divergence  at  40°  N,  consistent  with  the  climatological 
data.  March  through  June  (Figures  4.15c  through  f)  depict 
the  wind  divergence  moving  poleward  and  also  shows  the  wind 
speed  increasing  to  a  maximum  in  June,  especially  in  the 
southern  half  of  the  region.  Figures  4.15g  and  h  (July  and 
August)  depict  strong  winds  and  equatorward  flow  throughout 
the  model  domain.  By  September  (Figure  4.15i),  the  wind 
divergence  reappears  in  the  north  and  migrates  equatorward 
over  the  remaining  three  months  of  the  year  (Figures  4.15j 
to  1).  By  December  the  wind  divergence  is  at  37.5°  N 
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(Figure  4.151),  which  is  around  2.5°  farther  equatorward 
than  the  climatological  wind  divergence. 

As  in  Experiment  1,  a  relatively  weak  poleward 
coastal  surface  flow  ~2  cm  s'1  was  observed  in  the  poleward 
region  of  the  model  domain  (Figure  4.16).  In  response  to 
the  wind  pattern,  an  equatorward  surface  current  was 
observed  along  most  of  the  coast  by  model  day  90  (Figure 
4.17) .  A  weak  poleward  undercurrent  became  evident  by  model 
day  90  (Figure  4.18),  and  became  as  large  as  5  cm  s'1  in 
layer  5  (316  m)  by  day  210  (Figure  4.19). 

Upwelling  of  temperatures  cooler  than  15°  C  was 
observed  by  model  day  120  (as  observed  in  the  climatological 
run),  as  seen  in  Figure  4.20.  Meanders  in  the  surface 
current  also  formed  by  the  same  model  day  as  in  the 
climatological  run  (day  210)  as  depicted  in  the  velocity  and 
temperature  contour  plots  (Figures  4.21  and  4.22, 
respectively) .  By  model  day  270,  the  first  sign  of  eddies 
appear  (see  the  velocity  and  temperature  contour  plots  of 
Figures  4.23  and  4.24),  which  is  -30  days  later  than  in  the 
climatological  run.  The  amplitude  of  the  meander  of  the 
coastal  current  increases  to  the  south  as  a  result  of  the 
stronger  equatorward  winds  in  this  region.  Again,  these 
eddies  are  cold  core,  cyclonic,  and  about  30  km  in  diameter. 

Upwelling  of  the  coldest  surface  water  occurred  by 
model  day  240  (not  shown) ,  corresponding  to  the  period  of 
strongest  equatorward  winds.  Upwelling  did  not  disappear 
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until  model  day  330  (not  shown) ,  which  is  ~60  days  later 
than  the  climatological  results.  At  the  end  of  the  year, 
meanders  in  the  equatorward  surface  current  along  with 
embedded  eddies  remain  (Figure  4.25).  A  poleward 
undercurrent  is  also  discernible  at  the  end  of  the  year 
(Figure  4.26) .  A  very  weak  poleward  coastal  surface  current 
(not  shown)  has  reappeared  in  the  poleward  region  of  the 
model  domain.  As  in  Experiment  1,  no  filaments  developed 
during  this  experiment. 

Figure  4.27  shows  the  time  series  of  total  kinetic 
energy  per  unit  mass  over  the  entire  domain  for  1981  and 
shows  four  periods  when  the  kinetic  energy  becomes  quasi - 
steady,  model  days:  210  -  222,  240  -  252,  297  -  303,  and 
309  -  333.  Energy  transfer  plots  were  obtained  for  these 
periods  and  yielded  mixed  ( i . e . ,  both  barotropic  and 
baroclinic)  instabilities  in  all  cases  in  the  coastal, 
equatorward  region  of  the  model  domain.  These  results  are 
summarized  in  Table  4.1  (see  Appendix  B  for  the  energy 
transfer  plots  for  this  experiment) .  During  the  time  frame 
297  to  303,  barotropic  instability  (Figure  B.5)  was  much 
greater  than  baroclinic  (Figure  B.6)  instability. 

3.  Experiment  3  -  1983  Winds 

Experiment  3  represents  the  El  Niflo  year  of  1983  and 
was  forced  with  the  monthly  winds  depicted  in  Figure  4.28 
representing  January  through  December.  The  year  starts  off 
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with  predominantly  poleward  winds  for  the  first  three  months 
as  seen  in  Figures  4.28a  -  c.  April  shows  a  change  in  the 
wind  pattern  as  equatorward  winds  are  present  throughout 
(Figure  4.28d) .  These  winds  remain  equatorward  and 
strengthen  from  May  through  September  (Figures  4.28e  -  i) . 
The  winds  weaken  in  October  (Figure  4.28j) .  A  divergence  in 
the  wind  field  appears  near  38°  N  during  November  (Figure 
4.28k)  and  remains  through  December  (Figure  4.281) . 

In  response  to  the  three  consistent  months  of  strong 
poleward  winds  at  the  start  of  1983,  a  poleward  coastal 
surface  current  with  velocities  as  high  as  5  cm  s'1  is 
established  by  model  day  30  (Figure  4.29)  and  remains 
through  model  day  60  (not  shown) .  By  model  day  120  an 
equatorward  surface  current  develops  along  most  of  the  coast 
(Figure  4.30).  A  weak  poleward  undercurrent  is  present  at 
model  day  90  as  seen  in  Figure  4.31,  and  strengthens  to 
around  6  cm  s'1  by  model  day  195  (not  shown)  . 

Upwelling  and  meanders  both  appear  by  model  day  150 
(Figure  4.32).  The  meanders  are  also  discernible  in  the 
velocity  plot  (Figure  4.33).  Note  that  the  meanders  appear 
farther  poleward  than  in  the  previous  experiments .  Eddies 
start  to  form  by  model  day  240  as  seen  in  the  velocity 
(Figure  4.34)  and  temperature  (Figure  4.35)  fields.  As  in 
the  previous  experiments,  these  eddies  are  cold  core  and 
cyclonic. 


The  coldest  upwelling  temperatures  were  reached  by 
model  day  240  but  were  warmer  {only  14.0°  C  as  compared  to 
13.5°  C  and  11.5°  C  for  climatology  and  1981,  respectively) 
than  observed  in  previous  runs.  No  upwelling  was 
discernible  after  model  day  270.  By  the  end  of  the  year, 
meanders  in  the  equatorward  surface  current  along  with 
embedded  eddies  (Figure  4.36),  a  weak  poleward  coastal 
current  (not  shown) ,  and  a  poleward  undercurrent  (Figure 
4.37)  were  discernible.  As  in  the  previous  experiments ,  no 
filaments  developed  during  the  experiment. 

Figures  4.38a  and  b  show  the  time  series  of  total 
kinetic  energy  per  unit  mass  for  1983.  Quasi-steady  periods 
that  prompted  additional  evaluation  were,  model  days:  186  - 
192,  222  -  234,  249  -  258,  279  -  288,  312  -  315,  324  -  327, 
and  342  -  348.  The  results  for  the  instability  analysis  for 
this  year  are  summarized  in  Table  4.1  (see  Appendix  C  for 
the  energy  transfer  plots) .  Every  period  showed  mixed  (both 
barotropic  and  baroclinic)  instabilities  to  be  present.  The 
instability  analysis  for  model  day  279  to  288  was  the  only 
case  where  one  type  of  instability  was  noticeably  greater 
than  the  other,  and  in  this  case  barotropic  (Figure  C.7)  was 
greater  than  baroclinic  (Figure  C.8). 

4.  Summary  of  First  Year  Run  Results 

Table  4.2  summarizes  the  climatological,  1981,  and 
1983  experiment  runs  for  the  first  year.  In  all  first  year 
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runs,  both  a  surface  equatorward  current  and  a  poleward 
undercurrent  developed.  There  was  also  a  poleward  coastal 
surface  current  present  during  the  winter  season.  Upwelling 
was  a  seasonal  feature  which,  as  expected,  had  the  coolest 
temperatures  during  a  non-El  Nifio  year  (1981)  and  warmest 
temperatures  during  an  El  Nifio  year  (1983)  .  Both  meanders 
and  eddies  were  generated.  The  eddies  were  cold  core, 
cyclonic,  and  had  a  size  around  the  Rossby  radius  of 
deformation  (30  km) .  Along  with  an  equatorward  surface 
current,  both  meanders  and  eddies  remained  at  the  end  of 
each  year.  Instability  analyses  revealed  that  both 
barotropic  and  baroclinic  processes  were  important  eddy 
generation  mechanisms.  In  all  first  year  runs,  no  filaments 
were  observed. 

B.  EXPERIMENTS  SPUN  UP  FROM  THE  PREVIOUS  YEAR 
1.  Experiment  4  -  Climatological  Winds 

Experiment  4  was  forced  with  the  same  monthly 
averaged  climatological  winds  as  in  Experiment  1  (Figure 
4.1).  In  this  experiment,  the  model  was  continued  with  the 
spun-up  output  from  Experiment  1,  rather  them  starting  from 
rest.  For  convenience  this  experiment  is  referred  to  as 
*Climo  II"  because  climatological  winds  are  used  a  second 
time. 

In  response  to  the  wind  forcing,  an  equatorward 
surface  current  was  observed  by  model  day  396  (year  day  (yd) 
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30)  (Figure  4.39) .  Embedded  in  this  current  are  meanders 
and  cyclonic  eddies  (Figures  4.39  and  4.40) .  A  weak 
poleward  coastal  surface  current  was  also  present  initially 
(not  shown) .  A  poleward  undercurrent  was  established  by 
model  day  426  (yd  60)  (Figure  4.41) . 

Upwelling  was  observed  by  model  day  546  (yd  180) 
(Figure  4.42),  and  peaked  by  model  day  606  (yd  240)  (Figure 
4.43),  with  a  minimum  temperature  near  12°  C.  Cessation  of 
upwelling  occurred  by  model  day  726  (yd  360)  (not  shown) . 

In  previous  model  runs  anticyclonic  meanders  were 
observed,  but  they  did  not  close  off  to  form  anticyclonic 
eddies.  During  this  experiment,  in  addition  to  the 
predominantly  cyclonic  eddies,  a  warm  core  anticyclonic  eddy 
(defined  by  closed  contour  in  streamline  plot  -  not  shown) 
formed  by  model  day  456  (yd  90)  (Figure  4.44).  This 
anticyclonic  eddy  persisted  for  the  remainder  of  the  year. 

Figure  4.45  is  a  cross  section  in  the  middle  of  the 
domain  and  depicts  the  broad  meanders  in  the  CC  along  with  a 
strong  core  of  the  CUC  near  250  m  in  depth.  Between  model 
days  696  (yd  330)  and  726  (yd  360) ,  an  eddy  formed  as  the 
result  of  a  pinched-off  meander  (Figures  4.46  and  4.47).  A 
filament  was  observed  during  this  experiment  by  day  606  (yd 
240)  and  extended  about  300  km  offshore  (15°  C  isotherm) 
near  40°  N  (Figure  4.43)  . 

A  time  series  of  total  kinetic  energy  per  unit  mass 
over  the  entire  domain  is  shown  in  Figure  4.48.  Four 
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periods  were  investigated  further  because  the  kinetic  energy 
became  quasi-steady.  These  time  frames  correspond  to  model 
days:  378  -  408  (yd  12  -  42),  453  -  465  (yd  87  -  99),  504  - 

516  (yd  138  -  150),  and  624  -  636  (yd  258  -  270).  All  four 
of  these  periods  consisted  of  both  barotropic  and  baroclinic 
instabilities.  These  results  are  summarized  in  Table  4.3 
(see  Appendix  D  for  the  energy  transfer  plots  for  this 
experiment)  . 

2.  Experiment  5  -  1982  Winds 

Experiment  5  was  forced  with  monthly  winds  from 
1982.  This  experiment  continues  the  spun  up  case  using  the 
output  of  Experiment  2  (1981)  as  a  starting  point,  rather 
than  starting  from  rest.  Figure  4.49  shows  the  monthly  wind 
patterns  within  the  model  domain  starting  with  January  and 
running  through  December.  January  starts  with  the  typical 
divergent  flow,  this  year  slightly  to  the  north  at  42°  N 
(Figure  4.49a).  February  moves  the  divergence  equatorward 
to  40°  N  and  has  weaker  winds  to  the  south  (Figure  4.49b). 
The  flow  in  March  is  onshore  throughout  the  model  domain 
with  a  slight  poleward  component  in  the  north  (Figure 
4.49c).  April  has  the  divergence  at  42.5°  N  (Figure  4.49d). 
May  starts  the  all  equatorward  flow  and  is  the  strongest 
(Figure  4.49e),  this  regime  continues  in  June  (Figure  4.49f) 
through  October  (Figure  4.49i).  A  divergent  flow  returns  in 
October  at  42.5°  N  (Figure  4.49j),  and  migrates  equatorward 
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to  38.5°  N  in  November  (Figure  4.49k)  and  remains  in 
December  (Figure  4.491). 

As  in  the  previous  experiment,  the  year  starts  with 
a  fully  developed  equatorward  surface  current  which  contains 
embedded  meanders  and  eddies  (Figure  4.50).  Also  present  is 
an  anticyclonic  eddy,  and  a  poleward  surface  coastal  current 
in  the  northern  part  of  the  model  domain  (Figures  4.50  and 
4.51).  A  relatively  strong  (~9.5  cm  s*1)  poleward 
undercurrent  is  present  by  model  day  426  (yd  60) .  A  cross 
section  in  the  middle  of  the  domain  (Figure  4.52)  shows  this 
feature  along  with  a  strong  subsurface  flow  associated  with 
an  eddy  (Figure  4.53). 

Upwelling  began  by  model  day  516  (yd  150)  (Figure 
4.54),  and  peaked  by  model  day  606  (yd  240)  (Figure  4.55). 
The  coolest  water  observed  was  around  12.0°  C.  The 
upwelling  ceased  by  model  day  696  (yd  330)  (Figure  4.56). 

Four  cyclonic  eddies  and  one  anticyclonic  eddy  were 
presex-  (e.g..  Figure  4.53)  throughout  most  of  this 
experiment.  The  size  ranged  between  30  and  150  km  in 
diameter  and  the  westward  propagation  of  these  larger 
features  was  around  500  km.  As  in  Experiment  4,  a  filament 
was  observed  during  this  experiment.  By  day  600  (yd  234) 
(not  shown)  this  filament  extended  about  350  km  offshore 
near  38.5°  N.  A  cold  core  eddy  pinched  off  from  this 
filament  by  day  606  (yd  240)  (Figure  4.55). 
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A  time  series  of  total  kinetic  energy  per  unit  mass 
over  the  entire  domain  was  produced  (Figure  4.57),  and 
showed  two  quasi-steady  periods  for  additional  study.  These 
were  model  days:  507  -  516  (yd  141  -  150),  and  525  -  528 
(yd  159  -  162) .  As  in  the  previous  experiments,  the  energy 
transfer  plots  indicate  the  presence  of  both  barotropic  and 
baroclinic  instabilities.  Table  4.3  summarizes  these 
results  (see  Appendix  E  for  the  energy  transfer  plots  for 
this  experiment) . 

3.  Experiment  6  -  1984  Winds 

Experiment  6  was  forced  with  monthly  winds  from 
1984.  This  experiment  continues  the  spun-up  case  using  the 
output  of  Experiment  3  (1983)  as  a  starting  point,  rather 
than  starting  from  rest.  Figure  4.58  shows  the  monthly  wind 
patterns  within  the  model  domain  starting  with  January  and 
running  through  December.  January  has  a  divergent  flow  at 
43°  N  (Figure  4.58a) ,  while  February  shows  that  the 
divergence  has  moved  equatorward  to  39°  N  (Figure  4.58b). 
From  March  through  May  (Figures  4.58c  through  e) ,  the  wind 
divergence  moves  from  41°  N  to  43°  N  and  the  winds  become 
stronger.  June  is  the  start  of  equatorward  wind  flow 
(Figure  4.58f),  which  continues  through  July  and  August 
(Figures  4.58g  and  h) .  In  September,  the  wind  divergence 
reappears  in  the  north  at  47°  N  (Figure  4.58i) .  During 
October  and  November,  the  divergence  migrates  equatorward  to 
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37°  N  (Figures  4.58i  and  k) .  The  year  ends  in  December  with 
weak  winds  along  the  coast  with  a  divergence  at  40°  N,  and 
with  predominantly  equatorward  winds  further  off  the  coast 
(Figure  4.581) . 

At  the  start  of  the  experiment,  there  is  a  well- 
developed  equatorward  surface  flow  with  embedded  meanders 
and  eddies  (Figure  4.59).  There  is  also  a  weak  poleward 
coastal  surface  current  in  the  northern  corner  of  the  model 
domain  which  remains  unti3  model  day  426  (yd  60)  (Figure 
4.60).  A  relatively  strong  (~10  cm  s"1)  poleward 
undercurrent  is  evident,  as  seen  in  a  cross-section  in  the 
middle  of  the  model  domain,  at  model  day  546  (yd  180) 

(Figure  4.61) . 

Upwelling  begins  by  model  day  456  (yd  90)  (Figure 

4.62)  .  Corresponding  to  the  period  of  maximum  equatorward 
winds,  the  peak  occurs  at  model  day  576  (yd  210)  (Figure 

4.63)  in  which  cooler  temperatures  of  around  10°  C  are 
present  along  the  entire  coast.  This  year  had  upwelling  of 
cooler  water  earlier  than  previous  years  and  also  had  the 
coldest  water  which  was  a  result  of  the  strong  equatorward 
winds.  Cessation  of  the  cooler  water  occurs  by  model  day 
726  (yd  360)  (not  shown) . 

The  eddies  are  around  90  km  in  diameter  and  stay 
within  300  km  of  the  coast.  A  filament  occurs  by  model  day 
666  (yd  300)  (Figure  4.64).  This  filament  is  observed  near 
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38°  N  as  a  plume  of  15°  C  water  <*. id  extends  about  200  km  off 
the  coast. 

The  time  series  plot  of  total  kinetic  energy  per 
unit  mass  over  the  entire  domain  is  shown  in  Figure  4.65. 

The  kinetic  energy  is  still  fluctuating  because  the  model  is 
still  spinning  up  during  the  first  part  of  the  year.  Three 
periods  were  selected  for  additional  study,  model  days:  576 
-  591  (yd  210  -  225) ,  606  -  624  (yd  240  -  258) ,  and  645  - 
657  (yd  279  -  291) .  All  of  these  included  both  barotropic 
and  baroclinic  instabilities,  and  are  summarized  in  Table 
4.3  (see  Appendix  F  for  the  energy  transfer  plots  for  this 
experiment)  . 

4.  Summary  of  Second  Tear  Run  Results 

Table  4.4  summarizes  the  follow-on  (second  year) 
model  runs.  In  general,  a  weak  coastal  poleward  surface 
current  was  evident  during  the  winter  season  in  the  northern 
region  of  the  domain  for  all  years  (Climo  II,  1982  and 
1984) .  A  meandering  equatorward  surface  current,  already 
present,  at  the  start  of  these  experiments  was  maintained. 

A  poleward  undercurrent  was  generated  in  all  experiments. 

The  eddies,  already  present  at  the  start  of  the  experiments, 
were  both  cold  core,  cyclonic,  and  warm  core,  anticyclonic 
eddies  in  Climo  II  and  1982,  while  in  1984  there  were  only 
cold  core,  cyclonic  eddies.  The  size  and  offshore  extent  of 
the  eddies  were  consistent  with  the  first  experimental  runs. 
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Seasonal  upwelling  occurred  with  minimum  temperatures  of 
~12°  C  for  Climo  II  and  1982,  and  of  -10®  C  for  1984.  The 
year  1984  had  the  largest  extent  of  cool  upwelled  water  as 
expected  for  a  year  with  the  strongest  equatorward  winds. 
Filaments  developed  in  all  three  experiments,  but  was  most 
pronounced  (as  expected)  in  1984.  The  filaments  had  an 
offshore  extent  between  200  and  350  km.  Instability 
analyses  revealed  that  both  barotropic  and  baroclinic 
processes  were  involved  in  eddy  generation. 

5.  Experiment  7  -  Multiple  Climatological  Runs 

Experiment  7  was  an  extended  climatological  run  in 
which  the  climatological  winds  used  in  Experiment  1  (Figure 
4.1)  were  used  in  repetition  four  times.  The  third  and 
fourth  years  of  the  run  follow  Experiment  4  (Climo  II)  and 
are  dubbed  "Climo  III"  and  "Climo  IV",  respectively.  (The 
values  in  the  model  on  the  last  day  of  a  run  were  retained 
for  use  at  the  start  of  the  next  run,  rather  than  starting 
from  scratch  or  rest  as  is  done  in  a  cold  restart.) 

In  response  to  the  wind  forcing,  an  equatorward 
surface  current  was  observed  by  model  day  760  and  1125  (yd 
30)  for  Climo  III  and  Climo  IV,  respectively  (not  shown) . 
Embedded  in  this  current  were  meanders  and  cyclonic  eddies 
(not  shown) .  A  weak  poleward  coastal  surface  current  was 
also  present  initially  (not  shown) .  A  poleward  undercurrent 
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was  evident  early  on  by  model  day  820  and  1185  (yd  60)  for 
Climo  III  and  Climo  IV,  respectively  (not  shown) . 

Upwelling  was  observed  by  model  day  850  and  1215  (yd 
120)  for  Climo  III  and  Climo  rv,  respectively  (not  shown) . 
The  upwelling  peaked  by  model  day  970  (yd  240)  for  Climo  III 
(not  shown)  and  day  1305  (yd  210)  for  Climo  XV  (not  shown) , 
with  a  minimum  temperature  near  11°  C.  Cessation  of 
upwelling  occurred  by  model  day  1060  (yd  330)  for  Climo  III 
(not  shown) ,  and  was  still  present  in  small  patches  by  day 
1455  (yd  360)  for  Climo  IV  (not  shown) . 

Figure  4.66  and  4.67  show  the  total  kinetic  energy 
per  unit  mass  time  series  over  the  entire  domain  for  Climo 
III  and  Climo  IV,  respectively.  Climo  III  has  a  peak 
kinetic  energy  of  166  cm2s'2  around  model  day  915  (yd  185) 
while  Climo  rv  has  a  peak  of  145  cm2s‘2  at  model  day  1293 
(yd  198) .  Figure  4.68  depicts  the  total  kinetic  energy  per 
unit  mass  time  series  for  all  four  climatological  runs 
together.  It  clearly  shows  that  the  first  and  second  runs 
were  involved  in  the  generation  or  spin -up  of  features  in 
the  model,  while  the  third  and  fourth  runs  oscillated  about 
125  cm28~2  and  indicate  the  maintenance  of  these  features 
(i.e.,  a  quasi-steady  state). 

6.  Experiment  8  -  Multiple  Tears  (1981  -  1984) 

Experiment  8  was  a  run  of  two  additional  years  of 
the  model  using  the  monthly  winds  of  1983  (Figure  4.28)  and 
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1984  (Figure  4.58) .  This  experiment  differs  from  Experiment 
3  (1983)  which  was  started  from  scratch  (cold  restart) .  The 
runs  in  this  experiment  followed  in  succession  after 
Experiment  5  (1982) ,  and  used  the  values  in  the  model  on  the 
last  day  of  the  year  as  input  to  the  restart. 

In  response  to  the  wind  forcing,  an  equatorward 
surface  current  was  observed  by  model  day  760  and  1125  (yd 
30)  for  1983  and  1984,  respectively  (not  shown) .  Embedded 
in  this  current  were  meanders  and  cyclonic  eddies  (not 
shown) .  A  weak  poleward  coastal  surface  current  was  also 
present  initially  (not  shown) .  A  poleward  undercurrent  was 
evident  by  model  day  820  and  1185  (yd  60)  for  1983  and  1984, 
respectively  (not  shown) . 

Upwelling  was  observed  by  model  day  880  (yd  150)  for 
1983  (not  shown) ,  and  day  1155  (yd  60)  for  1984  (not  shown) . 
The  upwelling  peaked  by  model  day  940  and  1305  (yd  210)  for 
1983  and  1984,  respectively  (not  shown) .  The  minimum 
temperature  was  near  11®  C  for  both  1983  and  1984,  but  1983 
had  only  patches  of  this  water  while  1984  had  this  cool 
water  along  the  entire  coast.  Upwelling  was  still  present 
in  small  patches  by  day  1090  and  1455  (yd  360)  for  1983  and 
1984,  respectively  (not  shown). 

Figure  4.69  and  4.70  are  the  total  kinetic  energy 
per  unit  mass  time  series  over  the  entire  domain  for  1983 
and  1984,  respectively.  The  year  1983  has  a  peak  kinetic 
energy  of  128  cm2s'2  around  model  day  963  (yd  233)  and  1984 
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has  a  peak  of  215  cm2s'2  at  model  day  1368  (yd  273)  .  Figure 
4.71  depicts  the  runs  for  all  four  years  (1981  through  1984) 
together  on  one  total  kinetic  energy  per  unit  mass  time 
series.  It  appears  that  the  first  year  (1981)  as  seen  in 
Experiment  2  and  the  second  year  (1982)  as  seen  in 
Experiment  5,  were  involved  in  generation  or  spinning  up  of 
the  model.  The  third  year  (1983),  which  was  an  El  Nifio 
year,  had  an  energy  level  on  par  with  1982,  so  that  no 
increase  in  kinetic  energy  occurred.  In  the  fourth  year 
(1984) ,  the  energy  level  increase  significantly  to  the 
highest  value  (215  cm2s*2)  .  Comparing  this  figure  to  the 
climatological  time  series  (Figure  4.68)  for  the  last  two 
years  shows  that  the  maximum  (1984)  and  minimum  (1983) 
kinetic  energies  lie  on  either  side  (as  expected)  of  the 
climatological  year.  During  a  non-El  Nifio  year,  as  is  1984, 
features  such  as  equatorward  currents  and  upwelling  would 
strengthen  while  during  an  El  Nifio  year,  as  in  1983,  these 
features  would  be  weaker. 

7.  Summary  of  Multiple  Tear  Run  Results 

Table  4.5  summarizes  the  last  two  years  of  the 
extended  (four  year)  runs.  Listed  here  are  Climo  IV,  1983, 
and  1984.  Similar  results  are  obtained:  A  poleward  coastal 
surface  current,  a  surface  equatorward  current,  a  poleward 
undercurrent ,  along  with  embedded  meanders  and  eddies.  Of 
interest,  upwelling  occurs  later  in  1983  (at  yd  150)  than 
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for  the  other  years  (yd  60  to  yd  90) .  Although  the  minimum 
temperature  is  the  same  (11.0°  C)  for  all  years,  1984 
experienced  a  larger  area  of  cooler  water  (all  along  the 
coast)  than  that  of  the  El  Nifio  (1983)  year.  In  all 
experiments,  there  were  both  cold  and  warm  core  eddies. 
Finally,  filaments  were  observed  during  Climo  IV  and  1984, 
but  not  during  1983. 
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TABLE  4.1  SUMMARY  OF  INSTABILITIES  FOR  CLIMATOLOGICAL  YEAR 
(TOP) ,  1981  (CENTER) ,  AND  1983  (BOTTOM) . 


MODEL  DAYS 
for 

Climatology 

FEATURE  TYPE 
(Meander, 

Eddy,  Both) 

INSTABILITY 

(Barotropic, 

Baroclinic, 

Mixed) 

LOCATION 

(Coastal, 

Offshore, 

Both) 

207-216 

Meander 

Mixed 

Coastal 

228-231 

Meander 

Mixed 

Both 

297-303 

Both 

Mixed 

Both 

MODEL  DAYS 
for 

1981 

FEATURE  TYPE 
(Meander, 

Eddy,  Both) 

INSTABILITY 

(Barotropic, 

Baroclinic, 

Mixed) 

_____ 

LOCATION 

(Coastal, 

Offshore, 

Both) 

210-222 

Meander 

Mixed 

Coastal 

240-252 

Meander 

Mixed 

Both 

297-303 

Both 

Mixed 

Both 

309-333 

Both 

Mixed 

Both 

MODEL  DAYS 
for 

1983 

FEATURE  TYPE 
(Meander, 

Eddy,  Both) 

INSTABILITY 

(Barotropic, 

Baroclinic, 

Mixed) 

LOCATION 

(Coastal, 

Offshore, 

Both) 

186-192 

Meander 

Mixed 

Coastal  1 

222-234 

Meander 

Mixed 

Both  1 

249-258 

Meander 

Mixed 

Coastal 

279-288 

Meander 

Mixed 

Coastal 

312-315 

Both 

Mixed 

Both 

324-327 

Both 

Mixed 

Both 

342-348 

Both 

Mixed 

Both 
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TABLE  4.2  SUMMARY  OF  EXPERIMENTS  -  FIRST  YEAR. 


FEATURE 

CURRENTS! _ 

Sfc  Poleward 

Sfc  Bquatorward 

Undercurrent 
(poleward  only) 

UPWELLINGt 
(<15  °  C) _ 

Start _ 

Maximum 

(min  temp) _ 

Cessation _ 

MEANDERS! _ 

Start _ 

EDDIES ! 

Start _ 

Core _ 

Rotation _ 

Size  (diameter) 

FILAMENTS! 


330-360  days  330-360  days  330-360  days 

~90  days _ ~90  days _ ~120  days 

~60  days  ~120  days  -60  days 


-120  days _ -120  days _ -150  days 

-240  days  -240  days  -240  days 


(13.0°  C) _ (11.0°  C) 

-270  days _ -330  days 


(14.0°  C) 
-270  days 


-240  days _ |  -270  days 


cold 

cold 

cyclonic 

cyclonic 

-30-60  km 

-30-60  km 

-240  days 
cold 


Start 

Offshore 

extent 


n/a 


n/a 


n/a 


TABLE  4.3  SUMMARY  OP  INSTABILITIES  FOR  CLIMATOLOGICAL  II  YEAR 
(TOP) ,  1982  (CENTER) ,  AND  1984  (BOTTOM) . 


MODEL  DAYS 
for 

Climo  II 

FEATURE  TYPE 
(Meander, 

Eddy,  Both) 

INSTABILITY 

(Barotropic, 

Baroclinic, 

Mixed) 

LOCATION 

(Coastal, 

Offshore, 

Both) 

378-408 

(12-42) 

Both 

Mixed 

Both 

453-465 

(87-99) 

Meander 

Mixed 

Coastal 

504-516 

(138-150) 

Both 

wixed 

Both 

624-636 

(258-270) 

Both 

Mixed 

Both 

MODEL  DAYS 

FEATURE  TYPE 

INSTABILITY 

LOCATION 

for 

(Meander, 

(Barotropic, 

(Coastal, 

1982 

Eddy,  Both) 

Baroclinic, 

Offshore, 

Mixed) 

Both) 

507-516 

(141-150) 

Both 

Mixed 

Both 

525-528 

(159-162) 

Both 

Mixed 

Both  I 

MODEL  DAYS 
for 

1984 

FEATURE  TYPE 
(Meander, 

Eddy,  Both) 

INSTABILITY 

(Barotropic, 

Baroclinic, 

Mixed) 

LOCATION 

(Coastal, 

Offshore, 

Both) 

576-591 

(210-225) 

Meander 

Mixed 

Coastal 

606-624 

(240-258) 

Meander 

Mixed 

Coastal 

645-657 

(279-291) 

Meander 

Mixed 

Coastal 
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TABLE  4.4  SUMMARY  OF  EXPERIMENTS  -  SECOND  YEAR. 


t 


FEATURE 

CLIMO  II 

1982 

1984 

CURRENTS} 

Sfc  Poleward 

-366-426 
(1-60),  & 
696-726 
(330-360) 
days 

-366-426 
(1-60),  & 
696-726 
(330-360) 
days 

-366-426 
(1-60),  & 
696-726 
(330-360) 
days 

Sfc  Equatorward 

-396  days 
(30) 

-396  days 
(30) 

-396  days 
(30) 

Undercurrent 
(poleward  only) 

-426  days 
(60) 

-426  days 
(60) 

-546  days 
(180) 

UPWELLING: 

(<15  °  C) 

Start 

-546  days 
(180) 

-516  days 
(150) 

-456  days 
(90) 

Maximum 
(min  temp) 

-606  days 
(240) 

(12.0°  C) 

-606  days 
(240) 

(12.0°  C) 

-576  days 
(210) 

(10.0°  C) 

Cessation 

-726  days 
(360) 

-696  days 
(330) 

-726  days 
(360) 

MEANDERS: 

Start 

-396  days 
(30) 

-396  days 
(30) 

-396  days 
(30) 

EDDIES : 

Start 

-396  days 
(30) 

-396  days 
(30) 

-396  days 
(30) 

Core 

cold  <Sc  warm 

cold  &  warm 

cold 

Rotation 

cyclonic  & 
anticyclonic 

cyclonic  & 
anticyclonic 

cyclonic 

Size  (diameter) 

-30-90  km 

-30-150  km 

-30-90  km 

FILAMENTS: 

Start 

-606  days 
(240) 

-600  days 
(234) 

-666  days 
(300) 

Offshore 

extent 

-300  km 

-350  km 

-200  km 
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TABLE  4.5  SUMMARY  OF  EXPERIMENTS  -  MULTIPLE  YEARS. 


FEATURE 

1983 

1984 

CURRENTS: 

Sfc  Poleward 

-1096-1215 
(1-120),  & 
1395-1455 
(300-360) 
days 

-731-850 
(1-120),  & 
1060-1090 
(330-360) 
days 

-1096-1185 
(1-90),  & 
1425-1455 
(330-360) 
days 

Sfc  Equatorward 

-1125  days 
(30) 

-760  days 
(30) 

-1125  days 
(30) 

Undercurrent 
(poleward  only) 

-1155  days 
(60) 

-790  days 
(60) 

-1155  days 
(60) 

UP WELLING: 

(<15  0  C) 

Start 

-1185  days 
(90) 

-880  days 
(150) 

-1155  days 
(60) 

Maximum 
(min  temp) 

-1305  days 
(210) 

(11.0°  C) 

-940  days 
(210) 

(11.0°  C) 

-1305  days 
(210) 

(11.0°  C) 

Cessation 

-1455  days 
(360) 

-1090  days 
(360) 

-1455  days 
(360) 

MEANDERS: 

Start 

-1125  days 
(30) 

-760  days 
(30) 

-1125  days 
(30) 

EDDIES: 

Start 

-1125  days 
(30) 

-760  days 
(30) 

-1125  days 
(30) 

Core 

cold  &  warm 

cold  &  warm 

cold  &  warm 

Rotation 

cyclonic  & 
anticyclonic 

cyclonic  & 
anticyclonic 

cyclonic  & 
anticyclonic 

Size  (diameter) 

-30-80  km 

-30-80  km 

-30-80  km 

FILAMENTS: 

Start 

-1305  days 
(210) 

none 

-1305  days 
(210) 

Offshore 

extent 

-350  km 

n/a 

-350  km 
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Figure  4.1  Mid-monthly  averaged  climatological  (1980  - 
1989)  ECMWF  winds  in  m  s'1:  (a)  January,  (b)  February,  (c) 

March,  (d)  April,  (e)  May,  (f)  June,  (g)  July,  (h)  August, 
(i)  September,  (j)  October,  (k)  November,  and  (1)  December. 
The  latitudinal  grid  point  1  (129)  corresponds  to  35°  N 
(47.5°  N) ;  so  that  the  latitudinal  grid  point  of  52 
corresponds  to  40°  N.  Maximum  wind  vector  is  20  m  s'1. 
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Figure  4.1  (continued)  . 
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Figure  4.1  (continued) . 
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Figure  4.3  Experiment  1:  Surface  velocity  vectors  at  day 
90.  To  avoid  clutter,  velocity  vectors  are  plotted  at  every 
third  grid  point  in  both  the  cross- shore  and  alongshore 
direction,  and  velocities  less  than  5  cm  s'1  are  not 
plotted. 


Depth  (m) 


Distance  off  shore  (km) 


Figure  4.4  Experiment  1:  Cross-section  at  y-704  km 
(41.25°  N)  of  the  meridional  component  of  velocity  at  day 
60.  Dashed  lines  indicate  equatorward  flow,  while  solid 
lines  indicate  poleward  flow. 
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Figure  4.5  Experiment  1:  Cross-section  at  y~704  km 
(41.25°  N)  of  the  meridional  component  of  velocity  at  day 
90.  Dashed  lines  indicate  equatorward  flow,  while  solid 
lines  indicate  poleward  flow. 
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Temperature  contour  at  model  day  120.0 


Figure  4.6  Experiment  1:  Surface  temperature  at  day  120 
The  contour  interval  is  1°  C.  The  temperature  decreases 
towards  the  coast. 
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Figure  4.7  Experiment  1:  Surface  velocity  vectors  at  day 
210.  To  avoid  clutter,  velocity  vectors  are  plotted  at 
every  third  grid  point  in  both  the  cross -shore  and 
alongshore  direction,  and  velocities  less  than  5  cm  s'1  are 
not  plotted.  Note  start  of  meander  at  y~500  km. 
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Figure  4.8  Experiment  1:  Surface  temperature  at  day  210. 
The  contour  interval  is  1°  C.  The  temperature  decreases 
towards  the  coast.  Note  start  of  meander  at  y~500  km. 
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Temperature  contour  at  model  day  2-10.0 


Figure  4.9  Experiment  1:  Surface  temperature  at  day  240. 
The  contour  interval  is  1°  C.  The  temperature  decreases 
towards  the  coast.  Note  the  formation  of  cold  core  eddy  at 
y~230  km. 
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Figure  4.11  Experiment  1:  Surface  temperature  at  day  360. 
The  contour  interval  is  1°  C.  Note  meanders  and  eddies  are 
still  discernible  in  temperature  contours. 
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Figure  4.12  Experiment  1:  Surface  velocity  vectors  at  day 
360  as  in  Figure  4.3.  Note  equatorward  surface  current  over 
a  poleward  undercurrent.  The  number  associated  with  the  low 
at  y~800  km  corresponds  to  temperature ,  not  velocity. 
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Figure  4.13  Experiment  1:  Cross -sect ion  of  the  meridional 
component  of  velocity  at  day  360.  Note  meanders  and  eddies 
are  still  discernible  in  a  weak  equatorward  surface  current 
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Figure  4.14  Experiment  1:  Total  kinetic  energy  per  unit 
mass  time  series  for  climatological  (1980-1989)  year  over 
the  entire  domain  (Units  of  the  kinetic  energy  are  cm2s'2) 
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Figure  4.15  Mid-monthly  ECMWF  winds  for  1981  in  m  s'1: 

(a)  January,  (b)  February,  (c)  March,  (d)  April,  (e)  May, 
f  June,  (g)  July,  (h)  August,  (i)  September,  (j)  October, 
(k)  November,  and  (1)  December.  The  latitudinal  grid  point 
1  (129)  corresponds  to  35°  N  (47. 5°  N) ;  so  that  the 
latitudinal  grid  point  of  52  corresponds  to  40°  N.  Maximum 
wind  vector  is  20  m  s*1. 
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Figure  4.15 


(continued) . 
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Figure  4.15  (continued) . 


69 


Distance  off  shore  (km) 


1024.0  -853  3  -682  7  -512  0  -341.3 


1408.0 


1173.3 


234.7- 


Norlh-Soulh  ve!  contour  at  model  day  45.0 


Figure  4.16  Experiment  2:  Surface  meridional  component  of 
velocity  at  day  45.  The  contour  interval  is  1  cm  s'1.  The 
dashed  lines  indicate  equatorward  velocities.  The  maximum 
poleward  velocity  is  ~2  cm  s'1.  Latitudinal  alongshore 
distance  0  (1408)  km  corresponds  to  35°  N  (47.5°  N)  . 


Figure  4.17  Experiment  2:  Surface  velocity  vectors  at  day 
90.  To  avoid  clutter,  velocity  vectors  are  plotted  at  every 
third  grid  point  in  both  the  cross -shore  and  alongshore 
direction,  and  velocities  less  than  5  cm  s"1  are  not 
plotted. 
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Figure  4.18  Experiment  2:  Cress -section  at  y~704  km 
(41.25°  N)  of  the  meridional  component  of  velocity  at  day 
90.  Dashed  lines  indicate  equatorward  flow,  while  solid 
lines  indicate  poleward  flow.  Note  weak  poleward 
undercurrent . 
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Distance  alongshore  (km) 


Figure  4.19  Experiment  2:  Layer  5  (316  m)  velocity  vectors 
at  day  210.  The  maximum  poleward  velocity  is  5  cm  s'1  near 
y~469  km. 
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Temperature  contour  at  model  day  120.0 


Figure  4.20  Experiment  2:  Surface  temperature  at  day  120 
The  contour  interval  is  1°  c.  The  temperature  decreases 
towards  the  coast. 
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Figure  4.21  Experiment  2:  Surface  velocity  vectors  at  day 
210.  Note  the  meander  near  y-469  km. 
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Temperature  contour  at  model  day  210.0 


Figure  4.22  Experiment  2:  Surface  temperature  at  day  210 
The  contour  interval  is  1°  C.  The  temperature  decreases 
towards  the  coast.  Note  the  meander  near  y~469  km. 
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Figure  4.23  Experiment  2:  Surface  velocity  vectors  at  day 
270.  The  number  associated  with  the  high  at  y-150  km 
corresponds  to  temperature ,  not  velocity.  Note  the  eddy 
near  y~600  km. 
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Figure  4.24  Experiment  2:  Surface  temperature  at  day  270. 
The  contour  interval  is  1°  C.  The  temperature  decreases 
towards  the  coast.  Note  the  cold  core  eddy  near  y~6000  km. 
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Figure  4.25  Experiment  2:  Surface  temperature  at  day  360. 
The  contour  interval  is  1°  C.  The  temperature  decreases 
towards  the  coast.  Note  existence  of  meanders  and  embedded 
eddies . 
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Figure  4.26  Experiment  2:  Cross-section  at  y-704  km 
(41.25*  N)  of  the  meridional  component  of  velocity  at  day 
360.  Dashed  lines  indicate  equatorward  flow,  while  solid 
lines  indicate  poleward  flow.  Note  existence  of  poleward 
undercurrent . 
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Kinetic  Energy  vs.  Time 


Figure  4.27  Experiment  2:  Total  kinetic  energy  per  unit 
mass  time  series  for  1981  over  the  entire  domain  (Units  of 
the  kinetic  energy  are  cm2s'2) .  Quasi-steady  periods:  days 
210-222,  240-252,  297-303,  and  309-333. 
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Figure  4.28  Mid -monthly  ECMWF  winds  for  1983  in  m  s 
(a)  January,  (b)  February,  (c)  March,  (d)  April,  (e)  May, 

(f)  June,  (g)  July,  (h)  August,  (i)  September,  (])  October, 
(k.)  November,  and  (1)  December.  The  latitudinal  grid  point 
1  (129)  corresponds  to  35°  N  (47.5°  N) ;  so  that  the 
latitudinal  grid  point  of  52  corresponds  to  40°  N.  Maximum 
wind  vector  is  20  m  s*1. 
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Figure  4.28  (continued) . 
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Figure  4.28  (continued) . 
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Figure  4.30  Experiment  3:  Surface  velocity  vectors  at  day 
120.  To  avoid  clutter,  velocity  vectors  are  plotted  at 
every  third  grid  point  in  both  the  cross -shore  and 
alongshore  direction,  and  velocities  less  than  5  cm  s'1  are 
not  plotted.  Note  equatorward  surface  current  dominates. 
The  number  associated  with  the  high  at  y~250  km  corresponds 
to  temperature,  not  velocity. 
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Figure  4.31  Experiment  3:  Cross-section  at  y~704  km 
(41.25°  N)  of  the  meridional  component  of  velocity  at  day 
90.  Dashed  lines  indicate  equatorward  flow,  while  solid 
lines  indicate  poleward  flow. 
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Temperature  contour  at  model  day  150.0 


Figure  4.32  Experiment  3;  Surface  temperature  at  day  150 
The  contour  interval  is  1°  c.  The  temperature  decreases 
towards  the  coast.  Note  the  existence  of  upwelling  and 
meanders . 
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Figure  4.33  Experiment  3:  Surface  velocity  vectors  at  day 
150.  To  avoid  clutter,  velocity  vectors  are  plotted  at 
every  third  grid  point  in  both  the  cross -shore  and 
alongshore  direction,  and  velocities  less  than  5  cm  s'1  are 
not  plotted.  The  number  associated  with  the  high  at  y-850 
km  corresponds  to  temperature,  not  velocity.  Note  the 
meanders  near  y~938  km. 
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Figure  4.34  Experiment  3:  Surface  velocity  vectorB  at  day 
240.  To  avoid  clutter,  velocity  vectors  are  plotted  at 
every  third  grid  point  in  both  the  cross -shore  and 
alongshore  direction,  and  velocities  less  than  5  cm  s"1  are 
not  plotted.  The  number  associated  with  the  high  at  y~250 
km  corresponds  to  temperature,  not  velocity.  Note  existence 
of  eddies. 
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Figure  4.36  Experiment  3:  Surface  temperature  at  day  360. 
The  contour  interval  is  1°  C.  The  temperature  decreases 
towards  the  coast.  Note  eddies  remain  at  end  of  year. 
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Figure  4.37  Experiment  3:  Cross-section  at  y-704  km 
(41.25*  N)  of  the  meridional  component  of  velocity  at  day 
360.  Dashed  lines  indicate  equatorward  flow,  while  solid 
lines  indicate  poleward  flow.  Note  poleward  undercurrent 
remains  at  end  of  year. 
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Figure  4.38a  Experiment  3:  Total  kinetic  energy  per  unit 
mass  time  series  for  1983  (days  1-306)  over  the  entire 
domain  (Units  of  the  kinetic  energy  are  cm2s'2)  .  Quasi - 
steady  periods:  days  186-192,  222-234,  249-258,  and  279-288. 
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Figure  4.39  Experiment  4:  Surface  velocity  vectors  at  day 
396  (yd  30) .  To  avoid  clutter,  velocity  vectors  are  plotted 
at  every  third  grid  point  in  both  the  cross -shore  and 
alongshore  direction,  and  velocities  less  than  5  cm  s'1  are 
not  plotted.  The  number  associated  with  the  low  at  y~80  km 
corresponds  to  temperature,  not  velocity.  Note  equatorward 
surface  current  with  embedded  meanders  and  cyclonic  eddies. 
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Figure  4.40  Experiment  4:  Surface  temperature  at  day  396 
(yd  30) .  The  contour  interval  is  1®  C.  The  temperature 
decreases  towards  the  coast.  Note  cold  core  eddies. 
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Figure  4.41  Experiment  4:  Cross-section  at  y-704  km 
(41.25°  N)  of  the  meridional  component  of  velocity  at  day 
426  (yd  60) .  Dashed  lines  indicate  equatorward  flow,  while 
solid  lines  indicate  poleward  flow.  Note  existence  of 
poleward  undercurrent. 
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Temperature  contour  at  model  day  546.0 


Figure  4.42  Experiment  4:  Surface  temperature  at  day  546 
(yd  180) .  The  contour  interval  is  1°  C.  The  temperature 
decreases  towards  the  coast. 
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Figure  4.43  Experiment  4:  Surface  temperature  at  day  606 
(yd  240) .  The  contour  interval  is  1°  C.  The  temperature 
decreases  towards  the  coast. 
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Figure  4.44  Experiment  4:  Surface  velocity  vectors  at  day 
456  (yd  90) .  To  avoid  clutter,  velocity  vectors  are  plotted 
at  every  third  grid  point  in  both  the  cross -shore  and 
alongshore  direction,  and  velocities  less  than  5  cm  s'1  are 
not  plotted.  Note  anticyclonic  eddy  at  y~650  km. 
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Figure  4.45  Experiment  4:  Cross-section  at  y-704  km 
(41.25°  N)  of  the  meridional  component  of  velocity  at 
666  (yd  300) .  Dashed  lines  indicate  equatorward  flow 
solid  lines  indicate  poleward  flow. 
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Figure  4.46  Experiment  4:  Surface  temperature  at  day  696 
(yd  330) .  The  contour  interval  is  1°  C.  The  temperature 
decreases  towards  the  coast.  Note  formation  of  eddy  from 
pinched  off  meander  at  y~469  km. 
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Figure  4.48  Experiment  4:  Total  kinetic  energy  per  unit 
mass  time  series  for  Climo  II  over  the  entire  domain  (Units 
of  the  kinetic  energy  are  cm2s"2)  .  Quasi -steady  periods: 
days  378-408  (yd  12-42),  453-465  (yd  87-89),  504-516  (yd 
138-150),  and  624-636  (yd  258-270). 
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Figure  4.49  Mid-monthly  ECMWF  winds  for  1962  in  m  s'1: 

(a)  January,  (b)  February,  <c)  March,  (d)  April,  (e)  May, 

(f)  June,  (g)  July,  (h)  August,  (i)  September,  (j)  October, 
(k)  November,  and  (1)  December.  The  latitudinal  grid  point 
1  (129)  corresponds  to  35°  N  (47.5°  N) ;  so  that  the 
latitudinal  grid  point  of  52  corresponds  to  40°  N.  Maximum 
wind  vector  is  20  m  s*1. 
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Figure  4.49  (continued). 
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Figure  4.49  (continued) 
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Figure  4.50  Experiment  5:  Surface  velocity  vectors  at  day 
396  (yd  30) .  To  avoid  clutter,  velocity  vectors  are  plotted 
at  every  third  grid  point  in  both  the  cross -shore  and 
alongshore  direction,  and  velocities  less  than  5  cm  s'1  are 
not  plotted.  The  numbers  associated  with  the  lows  embedded 
in  the  meanders  and  eddies  along  the  coast  correspond  to 
temperature,  not  velocity.  Note  fully  developed  equatorward 
surface  current  with  embedded  meanders  and  eddies,  along 
with  poleward  surface  current  north  of  y~1173  km. 
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Figure  4.52  Experiment  5:  Cross-section  at  y~704  km 
(41.25®  N)  of  the  meridional  component  of  velocity  at  day 
426  (yd  60) .  Dashed  lines  indicate  equatorward  flow,  while 
solid  lines  indicate  poleward  flow.  Note  strong  poleward 
undercurrent  along  with  subsurface  flow  associated  with  an 
eddy. 
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Figure  4.53  Experiment  5:  Surface  velocity  vectors  at  day 
426  (yd  60) .  To  avoid  clutter,  velocity  vectors  are  plotted 
at  every  third  grid  point  in  both  the  cross -shore  and 
alongshore  direction,  and  velocities  less  than  5  cm  s'1  are 
not  plotted.  The  numbers  associated  with  the  lows  embedded 
in  the  meanders  and  eddies  along  the  coast  correspond  to 
temperature,  not  velocity. 
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Figure  4.54  Experiment  5:  Surface  temperature  at  day  516 
(yd  150) .  The  contour  interval  is  1°  C.  The  temperature 
decreases  towards  the  coast. 
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Temperature  contour  at  model  day  606.0 


Figure  4.55  Experiment  5:  Surface  temperature  at  day  606 
(yd  240)  .  The  contour  interval  if?  1°  C.  The  temperature 
decreases  towards  the  coast.  Note  maximum  extent  of 
upwelling. 
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Figure  4.58  Mid-monthly  ECMWF  winds  for  1984  in  m  s'1: 

(a)  January,  (b)  February,  (c)  March,  (d)  April,  (e)  May, 

(f)  June,  (g)  July,  (h)  August,  li)  September,  (j)  October, 

(k)  November,  and  (1)  December.  The  latitudinal  grid  point 
1  (129)  corresponds  to  35°  N  (47.5°  N) ;  so  that  the 
latitudinal  grid  point  of  52  corresponds  to  40°  N.  Maximum 
wind  vector  is  20  m  s'1. 
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Figure  4.58  (continued) . 
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Figure  4.59  Experiment  6:  Surface  velocity  vectors  at  day 
396  (yd  30) .  To  avoid  clutter,  velocity  vectors  are  plotted 
at  every  third  grid  point  in  both  the  cross -shore  and 
alongshore  direction,  and  velocities  less  than  5  cm  s'1  are 
not  plotted.  The  numbers  associated  with  the  low  and  high 
embedded  in  the  meanders  and  eddies  along  the  coast 
correspond  to  temperature,  not  velocity.  Note  equatorward 
surface  flow  with  embedded  meanders  and  eddies. 
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Distance  off  shore  (km) 
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North-South  vel  contour  at  model  day  426.0 


Figure  4.60  Experiment  6:  Surface  meridional  con^jonent  of 
velocity  at  day  426  (yd  60) .  The  contour  interval  is  5  cm 
s'  .  The  dashed  lines  indicate  equatorward  velocities. 
Latitudinal  alongshore  distance  0  (1408)  km  corresponds  to 
35°  N  (47.5°  N) .  Note  weak  poleward  flow  along  the  coast 
north  of  y~1273  km. 
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Figure  4.61  Experiment  6:  Cross-section  at  y-704  km 
(41.25°  N)  of  the  meridional  component  of  velocity  at  day 
546  (yd  180) .  Dashed  lines  indicate  equatorward  flow,  while 
solid  lines  indicate  poleward  flow.  Note  existence  of 
poleward  undercurrent. 
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Figure  4.62  Experiment  6:  Surface  temperature  at  day  456 
(yd  90) .  The  contour  interval  is  1°  C.  The  temperature 
decreases  towards  the  coast.  Note  start  of  upwelling. 


123 


124 


Distance  alongshore  (km) 


Distance  off  shore  (km) 
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Temperature  contour  at  model  day  668.0 


Figure  4.64  Experiment  6:  Surface  temperature  at  day  666 
(yd  300) .  The  contour  interval  is  1°  C.  The  temperature 
decreases  towards  the  coast.  Note  occurrence  of  filament  at 
y~250  km. 
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Figure  4.66  Experiment  7:  Total  kinetic  energy  per  unit 
mass  time  series  for  Climo  III  over  the  entire  domain  (Units 
of  the  kinetic  energy  are  cm2s*2)  .  Peak  value  of  166  cm2s'2 
at  day  915  (yd  185). 
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Figure  4.67  Experiment  7:  Total  kinetic  energy  per  unit 
mass  time  series  for  Climo  IV  over  the  entire  domain  (Units 
of  the  kinetic  energy  are  cm2s'2)  .  Peak  value  of  145  cm2s'2 
at  day  1293  (yd  198) . 
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Figure  4.68  Experiment  7:  Total  kinetic  energy  per  unit 
mass  time  series  for  Climo,  Climo  II,  Climo  III,  and  Climo 
IV  over  the  entire  domain  (Units  of  the  kinetic  energy  are 
cm2s'2)  . 
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Figure  4.71  Experiment  8:  Total  kinetic  energy  per  unit 
mass  time  series  for  1981,  1982,  1983,  and  1984  over  the 
entire  domain  (Units  of  the  kinetic  energy  are  cm2s'2)  . 
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V.  SUMMARY  AMD  RECOMMENDATIONS 


A.  SUMMARY 

The  different  experiment  runs  are  summarized  in  Tables 
4.2,  4.4,  and  4.5.  In  all  first  year  runs  (see  Table  4.2), 
both  a  surface  equatorward  current  and  a  poleward 
undercurrent  developed.  There  was  also  a  poleward  coastal 
surface  current  present  during  the  winter  season.  Upwelling 
was  a  seasonal  feature  which,  as  expected,  had  the  coolest 
temperatures  during  a  non- El  Nifio  year  (1981)  and  warmest 
temperatures  during  an  El  Niflo  year  (1983)  .  Both  meanders 
and  eddies  were  generated.  The  eddies  were  cold  core, 
cyclonic,  and  had  a  size  around  the  Rossby  radius  of 
deformation  (30  km) .  Along  with  an  equatorward  surface 
current,  both  meanders  and  eddies  remained  at  the  end  of 
each  year.  Instability  analyses  revealed  that  both 
barotropic  and  baroclinic  processes  were  important  eddy 
generation  mechanisms.  In  all  first  year  runs,  no  filaments 
were  observed. 

The  follow-on  (second  year)  model  runs  are  summarized  in 
Table  4.4.  In  general,  a  weak  coastal  poleward  surface 
current  was  evident  during  the  winter  season  in  the  northern 
region  of  the  domain  for  all  years  (Climo  II,  1982  and 
1984) .  A  meandering  equatorward  surface  current,  already 
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present,  at  the  start  of  these  experiments  was  maintained. 

A  poleward  undercurrent  was  generated  in  all  experiments. 

The  eddies,  already  present  at  the  start  of  the  experiments, 
were  both  cold  core,  cyclonic,  and  warm  core,  anticyclonic 
eddies  in  Climo  II  and  1982,  while  in  1984  there  were  only 
cold  core,  cyclonic  eddies.  The  size  and  offshore  extent  of 
the  eddies  were  consistent  with  the  first  experimental  runs. 
Seasonal  upwelling  occurred  with  minimum  temperatures  of 
~12°  C  for  Climo  II  and  1982,  and  of  ~10°  C  for  1984.  The 
year  1984  had  the  largest  extent  of  cool  upwelled  water  as 
expected  for  a  year  with  the  strongest  equatorward  winds. 
Filaments  developed  in  all  three  experiments,  but  was  most 
pronounced  (as  expected)  in  1984.  The  filaments  had  am 
offshore  extent  between  200  and  350  km.  Instability 
analyses  revealed  that  both  barotropic  amd  baroclinic 
processes  were  involved  in  eddy  generation. 

The  last  two  years  of  the  extended  (four  year)  runs  are 
summarized  in  Table  4.5.  Listed  here  are  Climo  IV,  1983, 
and  1984.  Similar  results  are  obtained:  A  poleward  coastal 
surface  current,  a  surface  equatorward  current,  a  poleward 
undercurrent,  along  with  embedded  meanders  amd  eddies.  Of 
interest,  upwelling  occurs  later  in  1983  (at  yd  150)  tham 
for  the  other  years  (yd  60  to  yd  90) .  Although  the  minimum 
temperature  is  the  same  (11.0°  C)  for  all  years,  1984 
experienced  a  larger  area  of  cooler  water  (all  along  the 
coast)  than  that  of  the  El  Nifio  (1983)  year.  In  all 
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experiments,  there  were  both  cold  and  warm  core  eddies. 
Finally,  filaments  were  observed  during  Climo  IV  and  1984, 
but  not  during  1983. 

Experiments  1  through  6  investigated  the  formation  or 
generation  of  features  such  as  currents,  meanders,  eddies, 
and  filaments.  In  Experiments  7  and  8  a  quasi- steady 
equilibrium  state  for  the  CCS  was  reached  and  the  features 
observed  to  develop  in  previous  years  were  maintained. 

Figure  5.1  plots  the  mean,  eddy,  and  total  kinetic 
energy  per  unit  mass  for  Climo,  1981,  and  1983  (Experiments 
1-3)  together.  During  these  first  generation  years. 

Figure  5.1a  shows  1981  had  the  most  total  kinetic  energy  and 
1983  (El  Nifio  year)  the  least.  This  also  holds  true  for 
both  the  mean  and  eddy  kinetic  energy  per  unit  mass  plots 
(Figures  5.1b  and  5.1c).  The  decrease  in  kinetic  energy 
(total  and  mean)  during  the  last  few  days  of  1983  is  not  an 
error,  but  rather  an  adjustment  for  continuation  with  the 
start  of  the  1984  kinetic  energy  plots. 

Figure  5.2  plots  the  mean,  eddy,  and  total  kinetic 
energy  per  unit  mass  for  Climo  II,  1982,  and  1984 
(Experiments  4-6)  together  (the  second  years  of 
generation).  The  largest  values  of  total,  mean,  and  eddy 
kinetic  energy  were  observed  during  1984  (Figure  5.2a 
through  5.2c) .  Climo  II  had  the  least  total  and  eddy 
kinetic  energy  (Figures  5.2a  and  5.2c),  but  both  Climo  II 
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and  1982  had  nearly  identical  low  mean  kinetic  energy 
(Figure  5.2b). 

In  all  experiments,  the  atmospheric  pressure  cycle  and 
subsequent  movement  of  low  and  high  pressure  patterns 
affected  the  wind  orientation  and  strength  over  the  model 
domain.  The  surface  current  was  observed  to  respond  to  the 
prevailing  wind  direction.  As  the  wind  turned  equatorward 
and  exerted  stress  on  the  water  surface,  Ekman  transport 
offshore  resulted  in  the  upwelling  of  cooler  water  along  the 
coast.  Interannual  variability  was  thus  observed  as 
expected. 

In  conclusion,  the  experiments  showed  that,  due  to  the 
wind  forcing,  features  such  as  currents,  meanders,  eddies 
and  filaments  could  be  generated  and  maintained.  These 
results  support  the  hypothesis  that  wind  forcing  is  a  very 
important  mechanism  for  the  generation  of  many  of  the 
observed  features  in  the  CCS. 

B.  COMPARISON  WITH  OBSERVATIONS  DURING  AN  EL  Nllfo  TEAR 

It  has  been  shown  that  the  model  results  qualitatively 
agree  with  available  observations  of  the  CCS  during  non-El 
NifLo  years.  For  example,  the  simulated  surface  and 
subsurface  currents,  along  with  eddies  and  filaments  agree 
in  size  and  magnitude  with  available  observations  of  the 
CCS,  as  listed  in  Batteen  at  al-  (1989) . 
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Here  the  focus  is  shifted  to  comparisons  of  observed  and 
model  results  during  an  El  Nifio  year.  Qualitative  results 
from  both  of  the  El  Nifio  experiments  (3  and  8)  along  with 
some  observations  taken  in  the  CCS  during  the  1983  El  Nifio 
are  presented  in  Table  5.1.  Because  the  model  does  not 
include  remote  forcing,  the  oceanic  teleconnection  will  not 
be  expressed  in  the  model  results.  There  is  still  good 
agreement  between  the  model  and  observations  taken  during 
the  1983  El  Nifio  within  the  model  domain.  In  particular, 
warmer  sea  surface  temperatures  along  with  onshore  advection 
are  evident,  especially  during  the  first  three  months  of  the 
year.  Poleward  surface  flow  (for  example  Figure  5.3)  is 
enhanced  and  the  entire  current  system  appears  weaker  them 
in  the  other  model  years.  Based  on  these  results,  it  is 
concluded  that  atmospheric  forcing  played  a  substantial  role 
in  the  anomalous  conditions  observed  in  the  CCS  during  the 
El  Nifio  of  1983. 

C.  NAVAL  RELEVANCE 

As  the  Navy  focus  has  recently  shifted  from  blue  water 
to  the  littoral,  EBC  models  such  as  the  one  used  here  take 
on  an  increased  importance .  Although  this  EBC  model  would 
require  a  finer  scale  before  it  could  be  used  to  predict  the 
location  of  features  of  interest,  it  is  a  start  in  the  right 
direction. 
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Accurately  predicting  the  location  of  meanders,  eddies 
and  regions  of  upwelling  becomes  critical  for  the  successful 
use  of  mines,  along  with  aviation,  amphibious,  and  covert 
swimmer  operations.  The  meanders  and  eddies  could  impact 
several  aspects  of  mining  (i.e.,  laying,  known  position,  and 
mine  removal) .  Smaller  amphibious  craft  could  get  swept  off 
course  if  caught  up  in  a  strong  jet  associated  with  a 
meander.  The  cold  upwelled  water  may  be  nutrient  rich  and 
support  a  plankton  bloom  which  could  significantly  reduce 
visibility  in  the  water.  In  addition,  the  colder  water 
adversely  affects  swimmer  immersion  times. 

Fog  often  results  in  upwelling  regions  where  wanner  air 
is  present  over  the  cooler  upwelled  water.  The  effect  of 
this  fog  is  to  reduce  surface  visibility,  which  could  be  a 
plus  (such  as  to  provide  cover  for  an  amphibious  landing)  or 
minus  (such  as  when  defending  an  area  against  covert 
operations) .  This  decreased  visibility  may  also  affect 
aviation  operations  by  requiring  instrument  flight  rules 
(IFR)  vice  visual  flight  rules  (VFR) .  Finally,  other 
optical  and  infrared  devices  will  be  degraded  by  the 
presence  of  increased  moisture  in  the  atmosphere  (i.e., 
fog)  . 

D.  RECOMMENDATIONS 

Future  studies  should  let  the  model  run  for  multiple 
years,  to  allow  for  spin-up.  As  seen  in  multiple  year 
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kinetic  energy  per  unit  mass  time  series  plots,  it  appears 
that  spin-up  may  take  only  three  to  four  years,  but  this 
needs  to  be  confirmed.  Additional  studies  should  use  daily 
rather  than  monthly  winds  to  examine  the  effects  of  events 
and  relaxations.  The  incorporation  of  bottom  topography  and 
a  more  realistic  coastline  would  allow  preferred  eddy 
generation  areas  to  be  identified.  The  incorporation  of 
bottom  topography  could  also  lead  to  an  even  better 
depiction  of  the  undercurrent,  which  plays  an  important  role 
in  eddy  generation.  Finally,  finer  horizontal  and  vertical 
resolution  should  be  incorporated  into  the  model  to  aid  in 
more  accurately  predicting  the  location  of  meanders,  eddies, 
and  filaments. 
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TABLE  5.1  COMPARISON  OP  OBSERVED  AND  MODEL  RESULTS  DURING 
AN  EL  NINO  YEAR  (1983) . 


El  Nifio  (1983)  observations 
in  CCS  compared  to  mean 

Experiment  3 
(1983  cold 
restart) 

Experiment  8 
(1983  extended 
run) 

Warmer  SST 
(1,  2,  3) 

Simulated 

Simulated 

Weaker  overall  current 
system 
(2) 

Simulated 

Simulated 

Enhanced  poleward  flow 
(2,  3) 

Simulated 

Simulated 

Enhanced  onshore  advection 
(1,  2,  3) 

Simulated 

Simulated 

Maximum  equatorward  surface 
current  farther  offshore 
(2) 

Simulated 

Simulated 

References:  (1)  Rienecker  and  Mooers  (1986) 

(2)  Simpson  (1983) 

(3)  Simpson  (1984) 


140 


KE  (cm  A2/sA2) 


Figure  5.1a  Experiments  1-3:  Total  kinetic  energy  per  unit 
mass  time  series  for  climatology  (1980-1989),  1981,  and  1983 
over  the  entire  domain  (Units  of  the  kinetic  energy  are 
cm2s*2)  . 
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Figure  5.1b  Experiments  1-3:  Mean  kinetic  energy  per  unit 
mass  time  series  for  climatology  (1980-1989),  1981,  and  1983 
over  the  entire  domain  (Units  of  the  kinetic  energy  are 
cm2s'2)  . 
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Figure  5.1c  Experiments  1-3:  Eddy  kinetic  energy  per  unit 
mass  time  series  for  climatology  (1980-1989),  1981,  and  1983 
over  the  entire  domain  (Units  of  the  kinetic  energy  are 
cm2s"2) . 
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Figure  5.2a  Experiments  4-6:  Total  kinetic  energy  per  unit 
mass  time  series  for  Climo  II,  1981,  and  1983  over  the 
entire  domain  (Units  of  the  kinetic  energy  are  cm2s'2)  . 
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Figure  5.2b  Experiments  4-6:  Mean  kinetic  energy  per  unit 
mass  time  series  for  Climo  II,  1981,  and  1983  over  the 
entire  domain  (Units  of  the  kinetic  energy  are  cm2s‘2)  . 
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Figure  5.2c  Experiments  4-6:  Eddy  kinetic  energy  per  unit 
mass  time  series  for  Climo  II,  1981,  and  1983  over  the 
entire  domain  (Units  of  the  kinetic  energy  are  cm2s'2)  . 
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Figure  5.3  Experiment  3:  Surface  meridional  component  of 
velocity  at  day  45.  The  contour  interval  is  1  cm  s'1.  Th« 
dashed  lines  indicate  equatorward  velocities.  Latitudinal 
alongshore  distance  0  (1408)  km  corresponds  to  35°  N  (47.5 
N) .  Note  enhanced  poleward  flow. 
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APPENDIX  A.  INSTABILITIES  FOR  EXPERIMENT  1  (CLIMO) 

This  appendix  contains  additional  analyses  from 
Experiment  1  (climatological  year)  for  the  periods  when  the 
total  kinetic  energy  per  unit  mass  became  quasi -steady. 
Energy  transfer  plots  which  consist  of  barotropic  (mean 
kinetic  energy  to  eddy  kinetic  energy)  and  baroclinic  (eddy 
potential  energy  to  eddy  kinetic  energy)  were  produced  for 
the  following  model  days:  207  -  216,  228  -  231,  and  297  - 
303.  These  plots  appear  as  Figures  A.l  through  A. 6. 
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Figure  A.l  Experiment  1:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for 
climatological  (1980-1989)  year,  for  model  days  207  to  216 
The  contour  interval  is  1  ergs  cm*3  s*1. 
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Figure  A. 2  Experiment  1:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  climatological  (1980-1989)  year,  for  model  days  207  to 
216.  The  contour  interval  is  1  ergs  cm'3  s'1. 
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Figure  A. 3  Experiment  1:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for 
climatological  (1980-1989)  year,  for  model  days  228  to  231. 
The  contour  interval  is  1  ergs  cm'3  s*1. 
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Figure  A. 4  Experiment  1:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  climatological  (1980-1989)  year,  for  model  days  228  to 
231.  The  contour  interval  is  l  ergs  cm'3  s'1. 
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Figure  A. 5  Experiment  1:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for 
climatological  (1980-1989)  year,  for  model  days  297  to  303. 
The  contour  interval  is  1  ergs  cm'3  s'1. 
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APPENDIX  B.  INSTABILITIES  FOR  EXPERIMENT  2  (1981) 


This  appendix  contains  additional  analyses  from 
Experiment  2  (1981)  for  the  periods  when  the  total  kinetic 
energy  per  unit  mass  became  quasi -steady.  Energy  transfer 
plots  which  consist  of  barotropic  (mean  kinetic  energy  to 
eddy  kinetic  energy)  and  baroclinic  (eddy  potential  energy 
to  eddy  kinetic  energy)  were  produced  for  the  following 
model  days:  210  -  222,  240  -  252,  297  -  303,  and  309  -  333. 
These  plots  appear  as  Figures  B.l  through  B.8. 
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Figure  B.l  Experiment  2:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.f  barotropic  energy  transfer)  for  1981, 
for  model  days  210  to  222.  The  contour  interval  is  10  ergs 


Distance  alongshore  (km) 


Figure  B.2  Experiment  2:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1981,  for  model  days  210  to  222.  The  contour  interval 
is  10  ergs  cm'3  s'1. 
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Figure  B.3  Experiment  2:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1981, 
for  model  days  240  to  252.  The  contour  interval  is  10  ergs 
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Figure  B.4  Experiment  2:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1981,  for  model  days  240  to  252.  The  contour  interval 
is  10  ergs  cm'3  s'1. 
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Figure  B.5  Experiment  2:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1981, 
for  model  days  297  to  303.  The  contour  interval  is  10  ergs 
cm'3  s'1. 
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Figure  B.6  Experiment  2:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1981,  for  model  days  297  to  303.  The  contour  interval 
is  10  ergs  cm'3  s'1. 
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Figure  B.7  Experiment  2:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1981, 
for  model  days  309  to  333.  The  contour  interval  is  10  ergs 
cm'3  s'1. 
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Figure  B.8  Experiment  2:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1981,  for  model  days  309  to  333.  The  contour  interval 
is  10  ergs  cm'3  s'1. 
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APPENDIX  C.  INSTABILITIES  FOR  EXPERIMENT  3  (1983) 


This  appendix  contains  additional  analyses  from 
Experiment  3  (1983)  for  the  periods  when  the  total  kinetic 
energy  per  unit  mass  became  quasi -steady.  Energy  transfer 
plots  which  consist  of  barotropic  (mean  kinetic  energy  to 
eddy  kinetic  energy)  and  baroclinic  (eddy  potential  energy 
to  eddy  kinetic  energy)  were  produced  for  the  following 
model  days:  186  -  192,  222  -  234,  249  -  258,  279  -  288,  312 
-  315,  324  -  327,  and  342  -  348.  These  plots  appear  as 
Figures  C.l  through  C.14. 
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Figure  C.l  Experiment  3:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1983, 
for  model  days  186  to  192.  The  contour  interval  is  10  ergs 
cm  3  s'1. 
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Figure  C.2  Experiment  3:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1983,  for  model  days  186  to  192.  The  contour  interval 
is  10  ergs  cm'3  s'1. 
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Figure  C.3  Experiment  3:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1983, 
for  model  days  222  to  234.  The  contour  interval  is  10  ergs 
-3  s'1 


Distance  off  shore  (km) 

18?4  8  853  3  -68? .  7  -512. §  -341  3  178  7  M 


Eddy  potential  energy  to  eddy  kinetic  energy 
average  over  model  days  222  0  to  234  0 


Figure  C.4  Experiment  3:  Energy  tranefere  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1963,  for  model  days  222  to  234.  The  contour  interval 
is  10  ergs  cm  *  s'1. 
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Figure  C.5  Experiment  3:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1983, 
for  model  days  249  to  258.  The  contour  interval  is  10  ergs 
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Figure  C.6  Experiment  3:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1983,  for  model  days  249  to  258.  The  contour  interval 
is  10  ergs  cm*3  s*1. 
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Figure  C.7  Experiment  3:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1983, 
for  model  days  279  to  288.  The  contour  interval  is  10  ergs 
cm'3  s'1. 
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Figure  C. 8  Experiment  3:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1983,  for  model  days  279  to  288.  The  contour  interval 
is  10  ergs  cm'3  s'1. 


175 


Distance  off  shore  (km) 

-1024.0  -853.3  -682.7  -512.8  -341.3  -170.7  0.0 


Mean  kinetic  energy  to  eddy  kinetic  energy 
average  over  model  days  312.0  to  315.0 


Figure  C.9  Experiment  3:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1983, 
for  model  days  312  to  315.  The  contour  interval  is  1  ergs 
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Figure  C.10  Experiment  3:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1983,  for  model  days  312  to  315.  The  contour  interval 
is  1  ergs  cm'3  s'1. 
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Figure  C.ll  Experiment  3:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1983, 
for  model  days  324  to  327.  The  contour  interval  is  1  ergs 
cm'3  s'1. 
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Figure  C.12  Experiment  3:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.f  baroclinic  energy  transfer) 
for  1983,  for  model  days  324  to  327.  The  contour  interval 
is  1  ergs  cm'3  s'1. 
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Figure  C.13  Experiment  3:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1983, 
for  model  days  342  to  348.  The  contour  interval  is  1  ergs 
cm  J  s  . 
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Figure  C.14  Experiment  3:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1983,  for  model  days  342  to  348.  The  contour  interval 
is  1  ergs  cm'3  s'1. 
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APPENDIX  D.  INSTABILITIES  FOR  EXPERIMENT  4  (CLHSO  II) 


This  appendix  contains  additional  analyses  from 
Experiment  4  (Climo  II)  for  the  periods  when  the  total 
kinetic  energy  per  unit  mass  became  quasi -steady.  Energy 
transfer  plots  which  consist  of  barotropic  (mean  kinetic 
energy  to  eddy  kinetic  energy)  and  baroclinic  (eddy 
potential  energy  to  eddy  kinetic  energy)  were  produced  for 
the  following  model  days:  378  -  408  (yd  12  -  42),  453  -  465 
(yd  87  -  99),  504  -  516  (yd  138  -  150),  and  624  -  636  (yd 
258  -  270).  These  plots  appear  as  Figures  D.l  through  D.8. 
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Figure  D.l  Experiment  4:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  Climo 
II,  for  model  days  378  to  408  (yd  12  -  42) .  The  contour 
interval  is  1  ergs  cm'3  s'1. 
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Figure  D.2  Experiment  4:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  Climo  II,  for  model  days  378  to  408  (yd  12  -  42) .  The 
contour  interval  is  1  ergs  cm'3  s'1. 
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Figure  D.3  Experiment  4:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  Climo 
II,  for  model  days  453  to  465  (yd  87  -  99) .  The  contour 
interval  is  1  ergs  cm'3  s'1. 
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Figure  D.4  Experiment  4:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  Climo  II,  for  model  days  453  to  465  (yd  87  -  99) .  The 
contour  interval  is  I  ergs  cm'3  s'1. 
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Figure  D.5  Experiment  4:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  Climo 
II,  for  model  days  504  to  516  (yd  138  -  150) .  The  contour 
interval  is  1  ergs  cm'3  s'1. 
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Figure  D.6  Experiment  4:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  Climo  II,  for  model  days  504  to  516  (yd  138  -  150) .  The 
contour  interval  is  1  ergs  cm'3  s*1. 
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Figure  D.7  Experiment  4:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  Climo 
II,  for  model  days  624  to  636  (yd  258  -  270) .  The  contour 
interval  is  1  ergs  cm"3  s"1. 
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Figure  0.8  Experiment  4:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  Climo  II,  for  model  days  624  to  636  (yd  258  -  270) .  The 
contour  interval  is  1  ergs  cm'3  s'1. 
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APPENDIX  B.  INSTABILITIES  FOR  EXPERIMENT  5  (1982) 


This  appendix  contains  additional  analyses  from 
Experiment  5  (1982)  for  the  periods  when  the  total  kinetic 
energy  per  unit  mass  became  quasi -steady.  Energy  transfer 
plots  which  consist  of  barotropic  (mean  kinetic  energy  to 
eddy  kinetic  energy)  and  baroclinic  (eddy  potential  energy 
to  eddy  kinetic  energy)  were  produced  for  the  following 
model  days:  507  -  516  (yd  141  -  150),  and  525  -  528  (yd  159 
-  162).  These  plots  appear  as  Figures  E.l  through  E.4. 
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Figure  E.l  Experiment  5:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.f  barotropic  energy  transfer)  for  1982, 
for  model  days  507  to  516  (yd  141  -  150) .  The  contour 
interval  is  1  ergs  cm"3  s'1. 
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Figure  E.2  Experiment  5:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1982,  for  model  days  507  to  516  (yd  141  -  150) .  The 
contour  interval  is  1  ergs  cm"3  s'1. 
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Figure  E.3  Experiment  5:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.f  barotropic  energy  transfer)  for  1982, 
for  model  days  525  to  528  (yd  159  -  162)  .  The  contour 
interval  is  1  ergs  cm'3  s'1. 
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Figure  E.4  Experiment  5:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1982,  for  model  days  525  to  528  (yd  159  -  162).  The 
contour  interval  is  1  ergs  cm'3  s*1. 
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APPENDIX  F.  INSTABILITIES  FOR  EXPERIMENT  6  (1984) 

This  appendix  contains  additional  analyses  from 
Experiment  6  (1984)  for  the  periods  when  the  total  kinetic 
energy  per  unit  mass  became  quasi -steady.  Energy  transfer 
plots  which  consist  of  barotropic  (mean  kinetic  energy  to 
eddy  kinetic  energy)  and  baroclinic  (eddy  potential  energy 
tc  eddy  kinetic  energy)  were  produced  for  the  following 
model  days:  576  -  591  (yd  210  -  225),  606  -  624  (yd  240  - 
258) ,  and  645  -  657  (yd  279  -  291) .  These  plots  appear  as 
Figures  F.l  through  F.6. 
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Figure  F.l  Experiment  6:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1984, 
for  model  days  576  to  591  (yd  210  -  225) .  The  contour 
interval  is  10  ergs  cm'3  s'  . 
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Figure  F.2  Experiment  6:  Energy  transfers  of  eddy  potential 
to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer) 
for  1984,  for  model  days  576  to  591  (yd  210  -  225) .  The 
contour  interval  is  10  ergs  cm"3  s'1. 
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Figure  F.3  Experiment  6:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1984, 
for  model  days  606  to  624  (yd  240  -  258) .  The  contour 
interval  is  10  ergs  cm*3  s*r. 
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Figure  F.5  Experiment  6:  Energy  transfers  of  mean  to  eddy 
kinetic  energy  (i.e.,  barotropic  energy  transfer)  for  1984, 
for  model  days  645  to  657  (yd  279  -  291) .  The  contour 
interval  is  10  ergs  cm'3  s'  . 


204 


Distance  alongshore  (km) 


206 


LIST  OF  RE7KRENCXS 


Allen,  J.S.,  1980:  Models  of  wind-driven  currents  on  the 
continental  shelf.  Ann.  .Rev.  Fluid  Mech.,  12,  389-433. 

Arakawa,  A.,  and  V.R.  Lamb,  1977:  Computational  design  of 
the  basic  dynamical  processes  of  the  UCLA  general 
circulation  model.  In,  Methods  in  Computational  Physics, 
J.  Chang,  ed..  Academic  Press,  17,  173-265. 

Batteen,  M.L.,  1989:  Model  simulations  of  a  coastal  jet  and 
undercurrent  in  the  presence  of  eddies  and  jets  in  the 
California  Current  System.  In,  Poleward  Flows  Along 
Eastern  Ocean  Boundaries,  S.J.  Neshyba,  C.N.K.  Mooers, 
R.L.  Smith  and  R.T.  Barber,  eds.,  Springer-Verlag,  263- 
279. 

Batteen,  M.L.,  and  Y.-J.  Han,  1981:  On  the  computational 

noise  of  finite-difference  schemes  used  in  ocean  models. 
Tellus,  33,  387-396. 

Batteen,  M.L.,  R.L.  Haney,  T.A.  Tielking,  and  P.G.  Renaud, 
1989:  A  numerical  study  of  wind  forcing  of  eddies  and 
jets  in  the  California  Current  System.  J.  Mar.  Res.,  47, 
493-523. 

Batteen,  M.L.,  M.J.  Rutherford,  and  E.J.  Bayler,  1992:  A 

numerical  study  of  wind-  and  thermal -forcing  effects  on 
the  ocean  circulation  off  Western  Australia.  J.  Phys. 
Oceanogr.,  22,  1406-1433. 

Bjerknes,  J.,  1966:  A  possible  response  of  the  atmospheric 
Hadley  circulation  to  equatorial  anomalies  of  ocean 
temperature.  Tellus,  18,  820-829. 

Blumberg,  A.F.,  and  G.L.  Mellor,  1987:  A  description  of  a 
three-dimensional  coastal  ocean  circulation  model.  In, 
Three-dimensional  Ocean  Models,  N.  Heaps,  ed.,  American 
Geophysical  Union,  4,  1-16. 

Breaker,  L.C.,  and  C.N.K.  Mooers,  1986:  Oceanic  variability 
off  the  central  California  coast.  Prog,  in  Ocesmogr., 

17,  61-135. 

Camerlengo,  A.L.,  and  J.J.  O'Brien,  1980:  Open  boundary 
conditions  in  rotating  fluids.  J.  Comput.  Physics,  35, 
12-35. 


207 


Carton,  J.A.,  1984:  Coastal  circulation  caused  by  an 
Isolated  storm.  J.  Phys.  Oceanogr . ,  14,  114-124. 

Carton,  J.A.,  and  S.6.H.  Philander,  1984.  Coastal  upwelling 
viewed  as  a  stochastic  phenomena.  J.  Phys.  Oceemogr., 

14,  1499-1509. 

Chelton,  D.B.,  1984:  Seasonal  variability  of  alongshore 

geostrophic  velocity  off  central  California.  J.  Geophys. 
Res.,  89,  3473-3486. 

Enfield,  D.B.,  1989:  El  Nifio,  past  and  present.  Rev. 
Geophys.,  27,  159-187. 

Haines,  R.T.,  1994:  A  numerical  study  of  interannual  wind 
forcing  effects  on  the  California  Current  System,  1980- 
1983.  M.S .  Thesis,  Naval  Postgraduate  School,  138  pp. 

Halliwell,  Jr.,  G.R.,  and  J.S.  Allen,  1987:  The  large-scale 
coastal  wind  field  along  the  west  coast  of  North 
America,  1981-1982.  J.  Geophys.  Res.,  92,  1861-1884. 

Han,  Y.-J.,  1975:  Numerical  simulation  of  mesoscale  eddies. 
PH.D.  Thesis,  University  of  California,  Los  Angeles,  154 
pp. 

Haney,  R.L.,  1974:  A  numerical  study  of  the  response  of  an 
idealized  ocean  to  large-scale  surface  heat  and  momentum 
flux.  J.  Phys.  Oceanogr.,  4,  145-167. 

Haney,  R.L.,  1985:  Midlatitude  sea  surface  temperature 

anomalies:  A  numerical  hindcast.  J.  Phys.  Oceanogr. ,  15, 
787-799. 

Hickey,  B.M.,  1979:  The  California  Current  System  - 
hypothesis  and  facts.  Prog.  Oceanogr.,  8,  191-279. 

Holland,  W.R.,  1978:  the  role  of  mesoscale  eddies  in  the 

general  circulation  of  the  ocean  -  numerical  experiments 
using  a  wind-driven  guasigeostrophic  model.  J.  Phys. 
Oceanogr.,  8,  363-392. 

Holland,  W.R.,  and  M.L.  Batteen,  1986:  The  parameterization 
of  subgrid  scale  heat  diffusion  in  eddy- resolved  ocean 
circulation  models.  J.  Phys.  Oceanogr.,  16,  200-206. 

Huyer,  A.,  1983:  Coastal  upwelling  in  the  California  Current 
System.  Prog,  in  Oceanogr.,  12,  259-284. 


208 


Huyer,  A.,  P.M.  Kosro,  S.J.  Lentz,  and  R.C.  Beardsley,  1989: 
Poleward  flow  in  the  California  Current  System.  In, 
Poleward  Flows  Along  Eastern  Ocean  Boundaries ,  S.J. 
Neshyba,  C.N.K.  Mooers,  R.L.  Smith  and  R.T.  Barber, 
eds..  Lecture  Notes  on  Coastal  and  Estuarine  Studies, 
Springer-Verlag,  142  - 159 . 

McCreary,  J.P.,  P.K.  Kundu,  and  S.-Y.  Chao,  1987:  On  the 

dynamics  of  the  California  Current  System.  J.  Mar.  Hes., 
45,  1-32. 

Mitchell,  R.P.,  1993:  A  numerical  study  of  seasonal  wind 
forcing  effects  on  the  California  Current  System.  M.S. 
Thesis,  Naval  Postgraduate  School,  125  pp. 

Mooers,  C.N.K. ,  and  A.R.  Robinson,  1984:  Turbulent  jets  and 
eddies  in  the  California  Current  and  inferred  cross - 
shore  transports.  Science,  223,  51-53. 

Nelson,  C.S.,  1977:  Wind  stress  and  wind  stress  curl  over 
the  California  Current.  NCAA  Tech.  Rep.  NMFS  SSFR-714, 
U.S.  Dept.  Commerce,  87  pp. 

O'Brien,  J.J.,  R.M.  Clancy,  A.J.  Clarke,  M.  Crepon,  R. 

Elsbery,  T.  Gammelsrod,  M.  MacVean,  L.P.  Roed,  and  J.D. 
Tompson,  1977:  Upwelling  in  the  ocean:  Two-  and  Three- 
dimensional  model  of  upper  ocean  dynamics  and 
variability.  In,  Modelling  and  Prediction  of  the  Upper 
Layers  of  the  Ocean,  E.B.  Kraus,  ed. ,  Pergamon  Press, 

New  York,  178-228. 

Pedlosky,  J.,  1974:  Longshore  currents,  upwelling  and  bottom 
topography.  J.  Phys.  Oceanogr.,  4,  214-226. 

Philander,  S.G.H.,  1990:  El  Ni&o,  La  Nlha,  and  the  Southern 
Oscillation.  Academic  Press,  Inc.,  289  pp. 

Philander,  S.G.H.,  and  J.-H.  Yoon,  1982:  Eastern  boundary 
currents  and  coastal  upwelling.  J.  Phys.  Oceanogr.,  12, 
862-879. 

Rienecker,  M.M. ,  and  C.N.K.  Mooers,  1986:  The  1982-1983  El 
Niflo  signal  off  northern  California.  J.  Geophys.  Res., 
91,  6597-6608. 

Semtner,  A. J. ,  and  Y.  Mintz,  1977:  Numerical  simulation  of 
the  Gulf  Stream  and  midocean  eddies.  J.  Phys.  Oceanogr., 
7,  208-230. 


209 


Simpson,  J.J.,  1983:  Large-scale  thermal  anomalies  in  the 
California  Current  during  the  1982-1983  El  Nifio. 

Geophys.  Res.  Lett.,  10,  937-940. 

Simpson,  J.J.,  1984:  El  Nifio- induced  onshore  transport  in 
the  California  Current  during  1982-1983.  Geophys.  Rea. 
Lett.,  11,  233-236. 

Tielking,  T.A.,  1988:  Wind  forcing  of  eddies  and  jets  in  the 
California  Current  System.  M.S.  Thesis,  Naval 
Postgraduate  School,  107  pp. 

Tisch,  T.D.,  S.R.  Ramp,  and  C.A.  Collins,  1992:  Observations 
of  the  geostrophic  current  and  water  mass 
characteristics  off  Point  Sur,  California,  from  May  1988 
through  November  1989.  J.  Geophya.  Rea.,  97,  12535- 
12555. 

Trenberth,  K.E.,  W.6.  Large,  and  J.6.  Olson,  1990:  The  mean 
annual  cycle  in  global  ocean  wind  stress.  J.  Phya. 
Oceanogr.,  20,  1742-1760. 

Weatherly,  G.L.,  1972:  A  study  of  the  bottom  boundary  layer 
of  the  Florida  Current.  J.  Phya.  Oceanogr.,  2,  54-72. 

Wickham,  J.B.,  A. A.  Bird,  and  C.N.K.  Mooers,  1987:  Mean  and 
variable  flow  over  the  central  California  continental 
Margin,  1978-1980.  Cont.  Shelf  Res.,  7,  827-849. 


210 


INITIAL  DISTRIBUTION  LIST 


No .  Copies 

1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  VA  22304-6145 

2.  Librarian,  Code  52  2 

Naval  Postgraduate  School 

Monterey,  CA  93943-5101 

3.  Chairman  (Code  OC/Bf)  1 

Department  of  Oceanography 

Naval  Postgraduate  School 
Monterey,  CA  93943-5122 

4.  Chairman  (Code  MR/Hy)  1 

Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey,  CA  93943-5114 

5.  Dr.  Mary  L.  Batteen,  (Code  OC/Bv)  2 

Department  of  Oceanography 

Naval  Postgraduate  School 
Monterey,  CA  93943-5122 

6.  Dr.  Curtis  A.  Collins,  (Code  OC/Co)  1 

Department  of  Oceanography 

Naval  Postgraduate  School 
Monterey,  CA  93943-5122 

7.  Dr.  Leslie  K.  Rosenfeld  (Code  OC/Ro)  1 

Department  of  Oceanography 

Naval  Postgraduate  School 
Monterey,  CA  93943-5122 

8.  LT  James  Vann  2 

12  Hedgerow  Place 

Jackson,  TN  38305 

9.  Director,  Naval  Oceanography  Division  1 

Naval  Observatory 

34th  and  Massachusetts  Ave,  NW 
Washington,  DC  20390 


211 


10 .  Commander  1 

Naval  Meteorology  and  Oceanography  Command 

1020  Balch  Boulevard 

Stennis  Space  Center,  MS  39529-5005 

11.  Commanding  Officer  1 

Naval  Oceanographic  Office 

Stennis  Space  Center,  MS  39529-5000 

12.  Commanding  Officer 

Fleet  Numerical  Meteorology  and  Oceanography  Center  1 
7  Grace  Hopper  Avenue,  STOP  1 
Monterey,  CA  93943-5001 


212 


